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CHAPTER I 
INTRODUCTION 


We have attempted to demonstrate the Wigner cusp a 
in neutron elastic scattering and total cross sections near the thresh- 
olds for excitation of the first excited states of the dominant iso- 
topes of Li, Fe, Zr, Ba, and Ce. The energy region investigated 
extended approximately 300 kev on either side of threshold. 


yt 


Theoretical work (see, for example, Newton Baz,~7 eee 
and pee? 1 provides detailed predictions of possible cusp effects. 
The effects originate essentially in two fundamental physical principles: 
the conservation of incident flux, expressed in the theory as unitarity 
of the scattering matrix; and the uncertainty principle of Heisenberg, 
expressed in the theory as analyticity of the scattering matrix. The 
magnitude and shape of cusp anomalies can yield information useful in 
scattering analysis and in establishing spins and parities of excited 
states. 

Cusp anomalies have been seen in proton elastic scattering near 


8,9, 10 


thresholds for (p,n) reactions. The theory shows that one essen- 


tial requirement for a sharp cusp is that coulomb effects be absent 
from that reaction channel which provides the onmeesioll Any reaction 


with outgoing neutrons, such as (p,n) (d,n) (a,n) (n,n'), satisfies 


this requirement. 





The advantage of studying Wigner cusps in elastic neutron, 
rather than elastic proton scattering, for example, is the easier 
analysis caused by the absence of coulomb effects in the incoming 
channel. After estimating possible cusp magnitudes, we concluded that 
the effects of inelastic neutron scattering reactions on elastic neu- 
tron scattering might be quite noticeable for selected cases. We 
hoped thereby to find a method for spin assignment of excited states 
and to elucidate interaction mechanisms between nuclei and neutrons 


of few Mev energy. 





CHAPTER II 
THEORY OF WIGNER CUSPS 


Previous articles on the theory of Wigner cusp phenomena have 
generally emphasized the simple spin cases (e.g., spinless particles); 
also, the results are not always given in a form readily compared with 
experiment. We would like to supplement the previous results by con- 
sidering more general spin and angular momentum configurations. Our 
results will provide a convenient form for experimental comparisons. 
In addition, the form of the results for pure Wigner cusp effects will 
be particularly suitable for subsequent deduction of "energy-averaged 


cusp expressions, as first derived in Ref. 6. 


Pure Wigner Cusp Theory 


We will compute the effect of reactions of the type X(n,n')X* 
on neutron differential elastic and total cross sections near threshold 
energies. Blatt and Pea canernie have derived necessary general ex- 
pressions for cross sections. For completeness, we give in Appendix A 
a summary of the physical basis for the cross section formulas, start- 
ing with the wave functions and the scattering matrix, 9, 

In applying these formulas we assume that elastic scattering 
with change of orbital angular momentum £ may be neglected. Because 
of the parity rule requiring Aé@=0, 2, 4... ., elastic scattering 


with change of £ turns out to be an unimportant phenomenon in most 
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Situations. In particular, in the present experiment this restric- 
tion is very mild, since, with one exception, the target nuclei were 
even-even. Hence, with I =0, i =3, s =% (I, i, and s are target, 
projectile, and channel spin, respectively), there is no possibility 
of elastic scattering with change of #£%. The exceptional case is that 
Or ti! (I = B/D) In this case, AL2Zec2 is limited by penetrability 
effects. For example, although incident #£ =0 waves theoretically 
can be scattered elastically into £#£ =2 outgoing waves (with com- 
pound spin J =1), it happens that in the energy region of interest 
kKR = 0.3 and any process involving £=2 is highly improbable. Con- 
sequently, for the present experiment, specializing the derivation to 
Ag =O is not a restriction on the validity of the results. 

If one considers angular momenta situations in which there can 
be scattering involving change of channel spin s, but no change in 
channel energy, complexities arise in a derivation of cusp effects. 

To preserve clarity in the presentation, we will restrict our deriva- 
tion to situations in which channel spin s is separately conserved, 
in addition to the previous assumption that orbital angular momentum 

£ is conserved. It will be shown that this restriction does not af- 
fect the validity of the derivation as it applies to present experi- 
mental results. We will indicate the nature of the additional complex- 
ities in cusp theory, which are present without this restriction, at 

a suitable point in the derivation. 

The following equations are numbered 3.16, 4.5, and 4.6 in 


Ref. le. We have specialized them to elastic scattering with no change 





of Be, 8, 6 ora. Becordingly, the dipremential elastic scattering 


eross section is given by 


Teta 
BO yg = = {(28 + 1)/[(2T + 1)(2i + v3} ee (1) 
s=|I-i| 
where 
Dy oo 
—4 2 e) 2 
donaweE LK /(2s + 1)] es B, (as ;a1s ) P, (cos yan (2) 
Ha@ecas) =(2/)) 5 dc 5 FC 
oe oe L 
i? 2 1 2 
x Z°( 2 Simic) (ey) 
ee eee 
a J 
al x 2 
x Re[(1-S.7, . )" (1-8 )] 
ash, ash, as £,, 308 £,, 


where dQ refers to the type of incident particle and the state of the 
struck nucleus; ae is en xX the deBrogtie reduced wave rengul or 
the incident neutron/target nucleon system; s, I, and i are channel, 
target, and projectile spins, respectively; and £, J are orbital and 
total angular momentum, respectively. Subscripts 1 and 2 are dumny 
indices for use in summation. All sums are unrestricted, except that, 


in practice, only one such sum, for example, J, runs to infinity, be- 


ale 
cause of the vanishing of %Z coefficients. Re stands for real part 
: J : 
of the bracketed expressior. "as £;s£ is the S-matrix element for 


elastic scattering with £ and s conserved (see Appendix A). 





L-£,+2£ 


1 3s 


Bs nl 
+ 1)4 (25, + 1)*(27, + 1) 


xs 
Qfat Ie (aie 0 


1) eet 4 


2(2. J, £, J 


spear ees 2 


ne 8 8 Jn» 8 L)(2, £ 


ial 2 0 Oj2, Z 


eo) 


2 zs 


1 


where W is the Racah eoefficient.!" Tables of 2 are available. 
The properties of W are discussed by Biedenharn, Blatt, and nose. 2” 

The energy dependence of the elastic cross section is contained 
in §. We will use the unitary and analytic properties of the S-matrix 
to obtain an expression for the elastic cross section in the vicinity 
of a reaction threshold. 

Among the inelastic reactions which can proceed from a given 
incident partial wave, we consider only the one with lowest threshold 
energy. Note, however, that an excitation of a higher excited state 
is not excluded as long as its threshold has the lowest energy of those 
reactions which may proceed from a particular incident partial wave. 


Baz? 


showed that the procedure for including higher inelastic thresh- 
olds of the same dominant partial wave is straight-forward in principle, 
although much more complicated. Since the present experiments were done 
in energy regions near first excited states, the case we consider is 
adequate, and the physical content of the results is not obscured by 
these complications. 

Also, other true inelastic processes (e.g., neutron capture) 
occurring at lower energy than the threshold energy of interest are 
neglected. 


In addition, we restrict the energy region in which the deriva- 


tion is valid to that in which the reaction is dominated by a single 


7? 





incident partial wave. In practice, such an energy region is usually 
sufficiently large to be of experimental interest. For example, an 
incident d wave contributes more than 90 percent of the cross section 
over an energy interval of 150 kev above threshold for inelastic neu- 
tron scattering leading to the 2" first excited state of wel One 
can explain this in terms of the optical model by noting that Ce has 

a large transmission coefficient for (low energy) outgoing Ss wave 
ie and that sufficiently close to threshold, inelastic neutron 
scattering is always dominated by that incident partial wave which is 
associated (through angular momentum conservation) with an outgoing 

Ss wave. 

Finally, to further simplify the presentation, but not as an 
essential limitation on the theory, we will only consider reactions 
dominated by outgoing s waves. The complete theory in this respect 
shows that reactions dominated near threshold by outgoing waves of or- 
bital angular momentum 4» greater than zero, produce cusp effects 
in the (2. )th derivative of the elastic cross seemen Therefore, 
these cases have less experimental interest. With one exception, this 
restriction is not essential since we measured threshold effects of 
reactions with outgoing s waves. The exceptional case is the reaction 
Ze Re ) ye (1.75 Mev) which Tucker, et al. have found dominated 


7 


by outgoing p waves already at about 30 kev above threshold. The 


90 


necessary modifications to the theory in the case of Zr” will be ob- 


vious. 


& = 





Under the above assumptions the reaction [see Eq. (5) below] 

may be described by its partial reaction cross section 
Rhee A (4 

where primes identify the particular channel quantities which couple 
to £', the dominant incident partial wave and a = (27' + 1)/(2s' + 1) 
is a statistical factor obtained after coupling angular momenta, averag- 
ing over initial magnetic quantum numbers and summing over final magnetic 
quantum numbers (see Blatt and eueslameiae or Blatt and pucaeaimi 

To simplify the notation, we now drop a number of subscripts. 
First @' is superfluous since we always consider a neutron incident 
upon the ground state of a target nucleus. Second, we reduce the double 
subscript set on the S-matrix element to a single set; the meaning of 
S, as given in Appendix A, should be kept in mind. 


Equation (4a) becomes: 


! t 


Ali 
L's 


) (4 ) 
Equation (4b) is essentially the statement that the scattering matrix 
S is univary (see Appendix A). The corresponding expression for the 


total reaction cross section is: 


. 1 ! 
a oe al L ie on 


Ce eee ee oT + ae Pres’ (5) 
S Eo i| Ze s'| 
Inverting Eq. (4b), the S-matrix element may be written: 


Tigi = (exp 24 87, ,)(2- 02), ,)? (6) 





1 
Equation (6) defines 285 10: as the (real) phase of the S- 
t 
matrix element Sins . To simplify the notation, we have expressed 
1 t 
oh in@units of xz - e . This normalization will be used in 


the remaining discussion except in instances where clarity requires 
that the units be specifically restored. 
Equation (6) is certainly valid above threshold. The usual 


2,3 


procedure is to expand the complete phase (which may have real and 


imaginary parts) of ee. 1 in energy about threshold. One consequence 


L's 
of this procedure is that the fixed value of the real phase at thresh- 
old appears in the final expressions. On the other hand, Eq. (6) is 
quite general and there is no reason to assume that the form of Eq. (6) 
should not persist even away from threshold (above threshold). Conse- 
quently, we shall allow the real part of the phase of S, i.e., 26, 

as well as O.> to vary over the energy region of interest, if the 
physical situation requires that they do so. This leaves us free later 
to examine the consequences of energy averaging over a cusp region where 
the phase 265 and the cross section on might be fluctuating rapidly 
because of the presence of many narrow resonances. 

Using R-matrix theory, heremnon! has shown that the S-matrix 
can always be written in a form equivalent to Eq. (6) if £, = 0, inde- 
pendently whether the energy is above or below threshold. An analytic 
continuation in energy across threshold always exists, even in the 
presence of many narrow and possibly overlapping resonances which would 
produce rapid variation with energy in the phase and in the reaction 


eross section. 
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The prescription for continuing Eq. (6) below threshold, assum- 
ing outgoing s wave neutrons in the reaction, is given by Kq. (29a) in 
Ref. 7; essentially, one may let —_ lo, above threshold go to 
ilo| below threshold. It is interesting that in R-matrix theory this 
prescription actually results from a change in the effective logarithmic 
derivative of the reaction channel wave function from a pure imaginary 
number above threshold to a pure real number below threshold. This is 
related to the fact that outgoing wave function (exp ik r)/r has to 
change to (exp -|k Ir) /r as the energy (E) is changed from atove 
threshold (E,,) to below. K | - JoM(E ie B,,) /a » M = reduced mass 
of the outgoing neutron. This prescription becomes more complicated if 
one considers outgoing waves in the reaction channel with ie 4 O. The 
necessary modifications for extension of the present theory to these 
cases is also given in Ref. 7. 


Using the fact that a” 


;_; << 1l near threshold, we can 
rs. 


simplify Eq. (6) to 


ae lees. - 
pig! ~ (exp Oi Sat) emcee) (7) 


and, upon substituting the analyticity prescription into Eq. (7), we 


obtain 


a | Sed iL ite 1 
Sptg! = (exp 21 gl = Fp gf) (8) 


' e 
where, following Newton, { 4} means A if EO? nan and B if 


EY = on where On is the reaction threshold energy. It should be 





dae 


J' ; ; J! 
noted that oO 51.1 in Eq. (8) is strictly equal to ote IE, BY) ; 


the accompanying {>} contains the analytic properties of the expres- 


eon . 
t 


:_+ into two terms. the first 
k's 


We have separated 5 


1 
exp ei St can be interpreted as the 'S-matrix element without 


g's! 3 
reaction effects, date L = o.! When this term is substituted into 
Eqs. (1), (2), and (3), and the indicated sums are performed, the re- 


sult is the "differential elastic cross section without reaction ef- 


feetey which we designate by dome - The second term, 
BA 
t t 1 
~(exp 2i By r1) 4 ae } produces a cusp term after substitution 


into the same equations. 


Accordingly, the Wigner cusp term kee is defined by the 


Pd 


set of equations 


do =do +Aaga (9) 
Cae a;a a;3a 
where 
T+i 
Ad = & oe i 
cara Ses 1)/(2) + 22 ee eee (10) 
Gite om a poe 
A dy6!.as' = [A /(2s' +1)] £ AB (as';as') P,(cos @)a2 (11) 


L=0 





A B, (as';as') Hone » : (hy J, a! ce ee) 


J , ' 
1 = ee, ee 


oe DY Z°(g' J’ 2 peels) 
rey, & — 
a 2 
J 
a ai) i ne 
x Re (om raat Brg) bog it (1 - exp ei ere) | J 
1 2 Ay fs 
where Opts! is still expressed in units of «tx G1 [see Eq. (eyi\., 


The sums over £_, J £ 


y Jy9 25» J. in Eq. (3) are reduced because of the 


a 


assumption that a single wave £' is dominant in the reaction. 


Because 2 £5 J, L, Jn. § L) = (aie 7. Lee Ji, 8 L) and 


se lee li 


Re(A* B) = Re(A B*), we may combine these sums, changing dummy indices 


a 


and dropping the extraneous index Q for simplicity, to obtain 


A B,(s';s') = (1/4) x : z 2 2° (4 Ne pte i Sly) 
a 


ee wang Tiger Tas ap 
x Re f(a = CAD 27. ye) (exp Qi S151) ) sir {| 


=(G ly) aes o°( 2 5 2 ee) 
a ae es 


— Je 1 ‘ai J | 
7 — 
x 2 SiH O95! 0 ye? Re aa expe L (285101 a 


Taking the Re [ ] in the last equation yields 
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A B,(s',s') = Sie S = z°( 4 TS ees" L) 
2 ie 

Jas J 

! sin(28 ’ ; eae 5 1 

x a eg ean a ; fs us (12) 

r£Zs Ls J J 

os(@6,,,_ = 5 9 

Ls Ls 


Within the limitations imposed during the derivation, Eqs. (9), 
(10), (11), and (12) describe expected Wigner cusp effects. The cusp 
effect is completely contained in the coefficients B, in Eq. (2); 


thee 1S, Bb. = Be 


, = Br +AB,, where AB is given by Eq. (12). 


Since it is frequent practice (cf. Ref. 20) to analyze neutron 
scattering data in terms of the Bo coefficients, the utility of Eq. 
(12) for cusp interpretation is clear. Should a cusp effect be found 
in only one or two of several experimental Be coefficients then the 
properties of the 24 coefficients severely restrict the participating 
partial waves. 

We now recall that the cusp predictions, Eas. (9) to (12) are 
not generally valid if scattering is possible which conserves energy 
but not channel spin, because most detection systems are only energy- 
selective. In such cases, the threshold reaction is coupled (via 
unitarity of the S-matrix) to other off-diagonal (different initial and 
final state channel spin) elements, as well as a diagonal element 
("pure" elastic scattering) of the S-matrix (see Eq. (4a)). An addi- 
tional sum over final state channel spins is required in Eqs. (1) to 


(3). As a result additional cusp terms appear in Eq. (11). Unless 


polarization is also measured, comparisons of theory and experiment 
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are ambiguous in such cases. None of the present experimental results 
required this extension to cusp theory for comparison purposes. 

The expected Wigner cusp in total cross sections (0,,) is 
easily obtained if we recall that, by assumption, the interaction in 


i 
theaeee)) D incident channel is purely elastic below threshold. 


n= J (ao/an)an + ona {ah 


Q 
I 


n° + E (A do/d2)an +o), tof | 


Mm 
Q 
+ 
> 
Q 


r 7 (13) 
where oe is the total cross section "without reaction effects." 
K 


Using the fact that [ P, (cos @)aN = 57 6 ist (where ss 
2 


Lae 


is a Kronecker delta) we obtain 


T+i = 
(25' +1 k 1 
= ee a) ha 


1 
Sabetituoine Zed £' J', 6' 0) = oF, Be i(-1)) > (2J +1)", where 
K K 
Oop! and Ory are Kronecker deltas, into the equation for AB. 


(with L=0), Eq. (12), and putting this result into Eq. (14), one 


obtains 
I+i Mig! 
Me = os en 1) a 7 meo(2T'+1) J! 
he aie | (a CLT (2s'+1) °re's' 
s'=|I-i| J'=|£'-s'"| 


- 2 she hae 7 
p } ‘ Seat toh (15) 


- sin eo 





1 


Oi gig: (see Eq. (6)). 


Finally, recognizing that the summations in Eqs. (15) and (5) are 


where we have explicitly indicated the units of 
identical, we may combine 


so 
{i}. {’ sin® 5» =" 
ie 
sin 28) tt 


Ss 


tll 
ae 
i Q 
O 
wn wn 
— 
af MN 
om 
NM & 
oO 
wm Cy 1 
Se 


and obtain the result 


° t 1 : 
PR a (25! + 1) - is x x (20 i) 
aa rarer) - + 1 r£'s' 
S mee i| a -S ‘| 
ey 
COs CB iat 
4 Et (16) 

- sin 25 p11 


Equation (16), combined with Eq. (13), is the general expression for 


Wigner cusp effect in the total cross section, within the limitations 


imposed by the derivation. 


Wigner Cusps--Spinless Particles 


We now extract the special case of spinless particles from these 


expressions, for comparison with other work. 


If all spins are zero, including that of the residual nuclei 


in the reaction, then Eq. (12) becomes 


- 0 sin(28, oy 
A B,(0;0) Sy aM ry? Meloy, © PR) 2 ae 5, my 


£ 0s(25, = ©» 
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el S 
Since Z(£ #200, O L) = (-1) (QL SReseO O OM) B:OomN0) 


K K 
i 1 
D or, » where © or, is a Kronecker delta, 


the sum in AB. is reduced to one term, and using Eqs. (10) and (11) 


L 


we obtain 


oo sin(28, - 5,) 
Mage ae 5 oy 2 sin By 4 2 be, an 
£=0 cos(25) - 8) 


£ 


where o) = 590° This expression can be written as 
il i * (ory xe) il 1 
A do = 2 Rel (5 1k : (1 - exp 2i be Ho dem 21 8, 
£=O0 
thus 
Ado = 2 Re |(1/.) £%(6) a exp 21 8 (o_) (17) 
2 O |i r’exp 


where f*(@) is the complex conjugate of the scattered amplitude, with 
its explicit form evident by comparison with the preceding line, and 
where (ae = (x — o5)3 i.e., the experimental reaction cross 
section. 

This is just the cusp term obtained by Bee for the case of 
spinless particles. Baz emphasized the explicit appearance of the am- 
plitude, suggesting that there is additional theoretical significance 
in cusp data beyond cross section measurements (squared amplitudes) 
that are usually made. We would like to caution that most cases in 


practice involve non-zero spins; consequently scattered amplitudes, 


in general, are azimuth angle (p) dependent. Comparison of theoretical 
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cusp effects with experiment may be undertaken if the theoretical ex- 

pressions have been averaged and summed over magnetic quantum numbers. 
This is essentially what we have accomplished, and, in general, little 
vesitge of complete amplitude forms is found in the resulting expres- 

sions. 


On the other hand, it is useful to see from Eq. (17) that one 


can estimate expected Wigner cusp magnitudes from the expression 
a 
JA do| = o (da)*/r 


the remainder of the cusp expression being of order unity. The ex- 
pected percentage cusp effect is inversely proportional to (ao)1/2 
and exploratory experiments might well be undertaken at angles where 
the cross section is least. 

The spinless particle result for total cross sections is a 
trivial reduction of Eq. (16) to 


cos 25, 
=g_° + 
On On (exp | ‘i 


- sin 255 


2 
This result is in agreement with those of others. ’ 
This concludes our discussion of pure Wigner cusp phenomena. 
Iwo other special cases, applicable to the present experiment, are de- 


duced in Appendix B. 
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Theory of "Energy-Averaged Wigner Cusps 


We turn now to the problem of calculating the expected threshold 
effect, if the measured cross sections are effectively energy averages 
of the actual detailed cross sections, containing closely spaced reson- 
ance structure. That the resulting effects are related intimately to 


Gh # 


pure Wigner cusp phenomena was first discussed by Meyerhof. 
showed that the energy average of Wigner cusp expressions leads to a 
predicted threshold effect which we shall call an "energy-averaged 
Wigner cusp. 

The physical situation we are considering is that the reaction 
and scattering cross sections are dominated by many closely spaced 
resonances. By closely spaced we mean a large number within the energy 
interval AE over which the averaging is performed: either experi- 
mentally, because of finite incident beam resolution or because of 
overlap of the resonances; or theoretically, by choosing AE>>D, 
where D is the average level spacing. No essential restric tenel 
will be placed on the level width to spacing ratio, I/D. 

Noting that the energy dependence of differential elastic cusp 


effects is contained in A Bes ‘;5") (see Eq. (12)) we can isolate the 


energy dependent quantities: 


t 
sin(28), ae 8 ) 
5 ed) 2 = 
( Opts! 2 sin ey 5 . by (18) 
cos(28 5101 ~ Date 


Brackets, ( - imply that an energy average is to be taken 


: see J 
over an interval containing many resonances. The phases Sos and 





1 


t 


ahi are rapidly varying with energy over the averaging interval and 
we will assume they are uncorrelated for J f a and/or £ # iplee MiGs 


is the well known statistical assumption for the compound nucleus. In- 
tuition might suggest that uncorrelated phases would lead immediately 
to breakup of expression (18) into a simple form, for example, the 


product of the averages of the factors. It turns out that such a break- 


up is by no means obvious, because of the presence of a within 
the averaging brackets. 

Using R-matrix theory and applying the random sign assumption 
to certain value quantities Yy» which will be discussed, Meyerhor! 


showed that 


and 
ais! 
- cos 28 
a oa 0 
( oi a i : ) a (20) 
mee sin os. p 
ts 
The value quantities vy), mentioned above were introduced by 
el 
Wigner and Kisenbud to parametrize the terms of the R-matrix, each 
of which has a resonance denominator. Thus, these quantities corre- 
spond to resonances and could be indexed in our notation by ae 
where 2X specifies a particular resonance of the set associated with 
angular momenta J, £ and spins. The random sign assumption leading 


to Eqs. (19) and (20) is made for value quantities of different index 
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% 45 of spe butmicenotejustifiedfor nelating valueequantities.of 
the game J, 2 and s; but different A; i.e., different resonances 
initiated by the same partial wave may have correlated phases. In 
particular, the reaction cross section appearing in the averaging 
brackets in Eqs. (19) and (20) proceeds through resonances correlated 
in phase to the remaining expression in the bracket, and it is surpris- 
ing that the simple results, Eqs. (19) and (20) were Ae Ce 

In Appendix C, we show how Eqs. (19) and (20) can be applied 
to expression (18) to give the result 
_ 1s" ) 
1 - 84.4) 


ri 
1 C 
ae a ate —_ 
Ls ! 
'e Ls ' 


One's! T 
0s(28) 


ae 
ms ris K K 
: } ee (21) 


where se is a Kronecker delta. 

Part of the proof of Eqs. (19) and (20), which led to Eq. (21), 
depends upon a cancellation of path integrals used in the energy averag- 
ing process, which will be complete only if normally slowly varying 
quantities remain slowly varying throughout the averaging interval. 

For the case of outgoing s waves in the reaction channel, one of these 
quantities (the logarithmic derivative of R-matrix theory) does not 
vary slowly right at reaction threshold. In essence, this difficulty 
is just the result of rediscovering pure Wigner cusp phenomena as one 


formally approaches very close to threshold in the averaging process. 
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This means that the derivation of Eqs. (19) and (20) is not strictly 
valid arbitrarily close to threshold. The actual experimental result 
would depend on the exact energy distribution of the beam. As long as 
the energy distribution of the beam does not overlap the threshold, 
the results, i.e., Eqs. (19) and (20) and consequently Eq. (21), will 
be applicable. 

After justifying Eq. (21) we have completed the derivation of 
"energy-averaged" cusp expressions. Recalling Eq. (12), which contains 
expression (18), and using Eq. (21), we obtain 


(go, t )- 
2 ré£s 
aig” _—- - iM t 1 t op s' 
( A B,(g 3s') ) (1/4) E. Lifes J L) 4 : y (22) 


the sums over £ and J have vanished because of Je 


in Eq. (21). The ‘energy-averaged'" differential cross section cusp term 


K 
and Oss! 


is obtained by substituting Eq. (22) in Eqs. (9-11). 
Substitution of relation (20) in Eq. (16) gives the expected 


result for the total cross section, 


(A o,) - { (25) 


That is, there is no "energy-averaged cusp in the total cross section. 
This is consistent with the basic optical model eeerecicas used in 
relating the model to experiment, that the energy-averaged total cross 
section, being linearly dependent upon the scattering matrix elements, 


is smooth in energy even across reaction thresholds. We will come back 


to this point when summarizing the theory. 
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"Rnergy-Averaged Cusps--Spinless Particles 


As an example, for comparisons to previous work, consider the 


special case of spinless bombarding and target particles 
i = 5 = 0 5 nde 


Equation (22) becomes 


( of ) 
a t t t t L'O 
( A B,(0;0) ere ewer" s' 2’ 2°50 L) | ; } (22a) 


2 2 2 
Weume Zee £' £'9£', OG) = (22' + 1g 2’ 0.0 | £' Bole, 
t 
indicating explicitly the units of Cory and by substituting Eq. 
(22a) into Eqs. (9-11), we obtain 


(Ado, =~ (22' +1)(x x-(24" + 1) (oF gs) in) me £' 0.0 eee L0)* 


yi Po dn 


- (28" + 1)(Co i) /de) (Ps) Lae (24) 


where Cie means units have been restored to (a ) and experi- 


re. 
mental values may be substituted for this quantity. Equation (24) is 


the same as the formula, Eq. (16) of Ref. 6. 
"Energy-Averaged" Cusps--Spin = Incident Particles 


We give the results for a second special case of Eq. (22) of 
particular interest in comparison with the present experimental data, 
where the incident particle is a neutron, the target nucleus is even- 


even, and incident d waves dominate the reaction. 
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Although not immediately obvious, Eq. (25) is identical to Eq. 


| Wscees = (£'+1)/(2'+1 


1 =2'xr/(e4''+1+ 2' 4r), where r is defined and used 
2 


(16') in Ref. 6. To compare, one must identify R 


+ g' Yr), Rot 


in Ref. 6. 
Since Eq. (25) is important to present experimental comparisons, 
it will be useful to give the terms in detail. It turns out that the 


cusp prediction, Eq. (25), is quite insensitive to the ratio Bion e 


recall that Ros was defined as the fraction of the reaction probabil- 
ity which originates from compound resonances of a particular J', for 


example £' +4. We make the reasonable choice R_, = 3(2J'+1)/(2e'+1) 


ar 
which is the same as the fraction of the total number of incident 


states of the age partial wave which have total angular momentum 


J'. With this choice of R Eq. (25) becomes 


J” 


(C0, 510) exp/ 043} 


= 
2 


Jie 
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[1 + 38/35 P,(cos @) + 18/35 P) (cos 6) Jan (26) 
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Equation (26) actually applies to any scattering experiments 
on medium and high atomic numbered even-even targets at energies where 
random phase concepts are useful, and near the threshold of a reaction 
dominated by incident d waves. The outgoing particle in the reaction 
channel must be a neutron, as mentioned previously, in order to obtain 
@ sharp cusp. 

In conclusion, we remark that the expressions for "energy- 
averaged” Wigner cusps depend only upon statistical assumptions, and 
not upon a particular model, for example, the optical model; rather, 
these expressions simply give the change in cross sections expected at 
reaction thresholds, due to unitarity of the scattering matrix. Meyerhof 
Bey however, that these statistical assumptions are consistent with 
the optical model assumption that the scattering matrix element consists 
of a smooth plus fluctuating Se 

Also, we note that "energy-averaged Wigner cusp expressions 
predict no below-threshold effect. Consequently, in regions of mass 
number and energy where "energy-averaged" theory applies, it is not 
necessary to specify the dominant outgoing wave in the reaction channel, 
in contrast to pure Wigner cusp theory (see the discussion just before 
Eq. (7)). In particular, outgoing p wave dominance of the reaction in 
the ee experiment, discussed earlier, causes no difficulty in analysis. 


Determination of the dominant incident partial wave in the reaction and 


Eq. (26) suffice. 


- 
—_ 





CHAPTER II1 
BXPERIMENT 


Method 


Neutron Differential Elastic Cross 
Section Measurements 


Neutrons were produced by the T(p,n) reaction, using tritium 
loaded zirconium and titanium targets for the proton beam from a High 
Voltage Engineering Company 3-Mev Van de Graaff generator stabilized 
to about 1 kev. Various tritium targets, described in Table I, of from 
10 to 70 kev thickness at 1.02 Mev proton energy (as determined by the 
rise method at threshold) were used. An eccentric circular motion was 
imposed upon the tritium target. Consequently, the proton beam struck 
a ring-shaped area on the target, which minimized cooling difficulties 
and also effectively averaged over nonuniform tritium deposition in 
the target, thereby producing a time integrated neutron flux very 
nearly proportional to accumulated proton charge. 

Tne neutron beam was collimated by a rectangular conical open- 


ing in LiCO.-loaded paraffin source shielding of characteristic thick- 


3 
ness 18 in. The geometry and shielding of the experiment is shown in 
Fig. 1. The collimation angle was such that the energy spread in the 
neutron beam due to T(p,n) kinematics was approximately 4 kev. A 
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standard long counter ~ placed at 0° with respect to the beam (see Fig. 


1) was used to monitor the neutron yield. 
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Fig. 1. Experimental Arrangement--Plan View: 1. Tritium 
loaded target; 2. Air-water spray coolant; 3. Eccentric wobble for 
additional cooling effect; 4. Paraffin shielding (loaded with 


Lic0,); 5. Paraffin (LiCO, loaded) rectangular conical collimator; 


5 
6. Scatterer in "In Beam" position; 7. Scatterer in "Out" position; 
8. Anthracene or stilbene crystal (a right circular cylinder, 2" 
dia x 1" long); 9. 6810A RCA photomultiplier tube; 10. Preampli- 


fier (linear and nonlinear gamma suppression outputs); ll. Teng 


Counter (BF) in monitor position. 
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The neutron beam was allowed to strike the scatterer, which was 
made slightly smaller than the aperture of the neutron beam, at about 
10 in. from the collimator end. The lithium scatterer was a right cir- 
cular cylinder oriented vertically; the remaining scatterers used were 
in the shape of a slab, oriented with their normals parallel to the 
floor and at 40° to the incident beam direction. The scatterers used 
are described in detail in Table II. Scatterer thicknesses were be- 
tween 1/5 and 5/8 neutron mean free paths at the energies of experi- 
mental interest. These rather thick scatterers were necessary to obtain 
sufficient counting rates. Scattered neutrons were detected at two 
angles by organic scintillation crystals (2 in. dia x 1 in. long anthra- 
cene and 2 in. dia x 3/4 in. long stilbene) placed one on each of the 
transmission and reflection sides of the scatterer at distances varying 
from 6 to 10 in. and with the crystals oriented as shown in Hage 1. 

The scatterer to crystal distance varied with detection angle and was 
set at the maximum consistent with a neutron packground/effect ratio 
less than two. A saturated space-charge pulse-shape discrimination 
output of the detection preamplifiers was used to suppress gamma-ray 
background. ‘Typical neutron pulse spectra, given in Fig. 2, illustrate 
that we also used an integral discriminator to bias out inelastic neu- 
trons. 

To facilitate neutron background subtraction, the scatterers 
were suspended on thin wires from an overhead crane so that the scat- 
terer could be alternately removed and replaced in the beam at each 


neutron energy (see Fig. 1). 





TABLE I 


Description of Tritium~-Loaded Targets 

















Index Made by Yield Thickness 
neutrons/sterad (at T(p,n) thresh- 
microcoulomb old 1.019 Mev) 
22050 Mev kev 
1 HAR 1.5 x 10° 22 
2 HAR O237 x 10° WA 
3 HAR ex 10° 60 
4 OR O33) x 10° DAE 
5 OR sale wax 10° aoa B 
HAR = Harwell, England (Tritium-loaded titanium) 
OR = Oak Ridge Natl. Lab., Tenn. (Tritium-loaded zirconium) 
TABLE II 
Description of Scatterers 
Scatterers Used in do Measurements 
Chemical Angle between 
Material Dimensions Purity slab and beam 
(inches ) (%) (degrees ) 
Cerium a. WO © 3-0 xX Ons 99 .8 50 
De 3.25 x 2.75 x 0.19 99.8 50 
Barium Pe OeexX 350. S< 219 97.5 hs 
Zirconium iO < 125 30-56 98 50 
Iron EO! ees ZS ee 00,52 98 50 
Lithium 3.05 dia X 4.4 high ~ 100 Rt. cir. cylinder 
Samples Used in 0,,( Transmission) Measurements 
Cerium iO) 6X 1.25 X 0.935 Nemec 
tf 8 ° 
Barium 1.0 dia. X 1.87 long coe gee, Cees 


sion was parallel 
tren 1.0 XxX 1.23) xX 0297 Hone to the beam axis. 











Fig. 2. Pulse spectra from detection equipment at neutron 
energy of 0.5 Mev. a. Gamma-Ray-Suppression Spectra (nonlinear pulses): 
1. With the detector (2" dia x 3/4" long stilbene crystal) in the neu- 
tron beam and with the pulse height analyzer (PHA) recording all preampli- 
fier nonlinear output pulses. 2. With the proton beam blocked by a beam 
stopper from hitting the tritium target, with a va source taped on the 
crystal, and with the PHA recording all preamplifier nonlinear output 
pulses. 3. With the proton beam on the tritium target, with the detector 
crystal in the neutron beam, and with a Na source on the crystal. Two 
integral discriminator biases were set, one on each of the linear and 
nonlinear spectra, and the PHA had been adjusted to accept nonlinear 
pulses which were associated with a time coincidence of discriminated 
linear and nonlinear pulses. b. Linear Spectra (linear preamplifier 
output pulses associated with a time coincidence of discriminated linear 
and nonlinear pulses): A. With the detector crystal in a typical experi- 
mental position (90° lab) and the scatterer in the experimental position 
in the neutron beam. B. With the detector in the same experimental posi- 
tion as in A, and with the scatterer retracted to its "Out" position 


(see Fig. 1). C. (A - B) = elastically scattered neutron pulse spectrun. 
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Neutron Total Cross Section Measurements 


The same neutron beam, source shielding, collimation, and gamma 
suppression described above were used in measuring total cross sections. 

Neutrons were detected in this case by a stilbene crystal with 
dimensions 1/2 in. dia X 1/2 in. long mounted at 0° to the incident 
beam. Transmission samples were right parallelepipeds whose area normal 
to the beam (* 1 sq. in.) was sufficient to shield the crystal from 
direct neutrons. (Neutrons were produced over a proton beam spot area 
of about 1/4 NE gt Pd 1/4 in.) The samples were 10 in. from the collima- 
tion (28 in. from the neutron source) and 20 in. from the crystal. 
Sample thicknesses were chosen to allow approximately 60 percent neutron 
transmission at energies of interest in the aepemmene.°" The sample 
and a Qin. long 1 in. dia. tungsten shield were mounted on thin sup- 
porting rods above the ends of two arms of a rotatable horizontal cross. 
A third arm had only a thin supporting rod attached. Thus it was pos- 
sible to count transmitted, background, and direct neutrons in sequence 


at each neutron energy by rotating the cross. 


Apparatus 


It remains to describe the electronics systems necessary to re- 
duce crystal scintillation events to (elastically scattered) neutron 
counts. A detailed discussion of our gamma suppression circuit is con- 
tained in Appendix D. The block diagram, Fig. 3, will be useful in 


following the present discussion. 
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Fig. 3. Block diagram of electronics, showing the function of 


each major component. 
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For each of two detection systems an RCA 6810A photomultiplier 
and preamplifier (Fig. 22) produced linear and nonlinear (gamma sup- 
pression) outputs in response to a crystal scintillation. These two 
outputs fed into a double delay line and RC coupled linear amplifier, 
respectively. 

Each of the amplifier outputs drove an adjustable bias circuit 
(integral discriminator) which generated a standard square 20 volt 
output for each input pulse that exceeded the bias level. After time 
synchronization was established by the use of a suitable delay cable, 
these pulses were allowed to enter a standard coincidence unit. A co- 
incidence output pulse was obtained for each simultaneous occurrence 
of input pulses. Thus, with properly adjusted bias levels on each of 
the linear and nonlinear pulses, a coincidence output pulse signalled 
the detection of an elastically scattered (linear pulse above linear 
bias) neutron (nonlinear pulse above nonlinear bias). These pulses 
were the experimental counts. 

A 256-channel pulse-height analyzer (PHA) was a valuable tool 
in setting up the gamma suppression circuit, the amplifier gains, and 
the discriminator biases. We also used the PHA as a monitor on the 
overall electronics performance during a run; these monitor spectra had 
a specific role in data analysis for some of the experiments (to be dis- 
cussed). Typical pulse spectra are shown in Fig. 2. The coincidence 
output pulses also served to gate the PHA when spectra such as those 
shown in Fig. 2 were taken. 


2 
A standard Hansen and McKibben long counter 3 was employed to 
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monitor the neutron yield. Direct neutrons were moderated by a paraffin 
cylinder and some were captured by wou in a one-in. dia. BF detec- 
tion tube. The subsequently emitted alpha particle initiated a pulse 
which was amplified and counted as a neutron event. The efficiency of 
the counter was approximately 2 xX 10°? and was fairly constant for a 
range of neutron energies from 0.5 Mev to 3 Mev. A 1/8 in. cadmium and 
3 in. paraffin shield on the sides and back of the counter reduced 

floor scattered neutrons to less than 15 percent of direct neutrons. 
Although absolute neutron yields of our tritium targets were measured 
using the long counter (see Table I), these yields were not required 


in the data analysis. 
Procedure 


Differential Elastic Cross Section 


Since the experimental interest centered on effects near the 
energies of first excited states, initial adjustments of the apparatus 
were made at these energies. For each experiment the nonlinear preamp- 
lifier controls (see Appendix D) were adjusted for optimum gamma-ray 
suppression, while the crystal was located in the scatterer "in beam" 
position, later called the "calibration position." Care was taken to 
insure that the crystal plane was parallel to the incident beam axis. 
This aspect corresponded to the aspect of the crystal toward scattered 
neutrons when it was in experimental position (see Fig. 1). After 
these adjustments were completed, the neutron yield and overall effi- 


ciency of the erystal/electronics system were calibrated by counting 
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direct target neutrons at several neutron energies above and below 
threshold. This was done for each of two crystal detector systems. 

After the crystals were moved to the experimental position 
(Fig. 1), the same energy range was covered in energy increments of 
approximately cne-half incident neutron-energy spread. At each energy, 
counts were recorded both with the scatterer in and with the scatterer 
out of the beam, thus providing a basis for subtraction of neutron 
background. On a succeeding run down in energy the same increments 
were used, but the points interleaved those from the run up; thus, we 
checked data reproducibility. 

All counts were normalized to a fixed amount of accumulated 
proton charge. When using 30-50 microamperes proton current, we found 
that the necessary counting time was typically 10 minutes for each 
energy point. 

Long counter monitor neutron. counts were also recorded at each 
point. Monitor neutron counts with the scatterer out of the beam 
served both as a check on the condition of the tritium target, which 
was subjected at times to as much as 70 watts of heat flow, and as the 
basis for a small correction of the data analysis for target yield 


Variation. 


Total Cross Sections 


As stated above, the gamma suppression circuit was adjusted at 
a neutron energy corresponding to an excited state of interest. After 


carefully aligning the sample and shield between the proton beam spot 





38 


and crystal detector, we took a sequence of counts at each energy point. 
By rotating the supporting horizontal cross on which sample and shield 
were mounted we were able to place first the sample, then the shield, 
and last the duplicate sample support (for a count of direct neutrons) 
into the beam. After each such set of counts the neutron beam energy 
was increased by approximately one-half incident beam-energy spread 
until a desired energy range was covered below and above threshold. A 
subsequent run down in energy, with interleaving points, verified re- 


producibility. 





CHAPTER IV 


DATA REDUCTION AND ANALYSIS 


Data Normalization 


The basic normalization used in the experiment was to a fixed 
accumulated proton charge, which implies an assumption that the result- 
ing time-integrated neutron flux was also a fixed quantity, except for 
the smooth change in the T(p,n) reaction cross section with energy, 
and except for a +1 percent uncertainty introduced by our proton cur- 
rent integrator. As mentioned before, the eccentric circular motion 
of the tritium target minimized any imperfection in this assumption, 
by averaging over nonuniformities of tritium deposition. However, 
small variations in integrated neutron flux from the true smooth shape 
of the T(p,n) reaction cross section might have been present when the 
neutron energy was changed, because of changes in proton beam intensity 
distribution which occur when one sets the accelerator to a new energy. 
If corrections were not made this could have produced scatter in the 
data, as plotted versus energy, in excess of statistical fluctuations. 
For this reason, monitor neutron counts, which had very small statis- 
tical errors, were used to renormalize the differential cross section 
effect counts to a smooth curve drawn through monitor neutron data; 
this curve was assumed to have the true smooth shape of the neutron 


yield versus neutron energy. The normalization was then a percentage 


Sy) 
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correction to each data point depending upon the percent deviation of 
the monitor counts from the smooth curve. This correction was usually 
less than one percent and never more than three percent. It was gen- 
erally impossible to distinguish by eye the difference between corrected 
and uncorrected data, but in seeking cusp effects, we wished to take 

all precautions necessary to establish the relative magnitude of the 


data points plotted versus neutron energy. 


Data Averaging 


The original data for experiments on zirconium, cerium and 
barium showed reproducible fluctuations in addition to a well defined 
drop in magnitude over a broad region above reaction thresholds. It 
turned out that energy averaging the original data did not mask the 
essential structure or features. Figure 4 illustrates this point, 
where we have used differential elastic cross section data from a 
cerium experiment. Therefore, we were able to improve the statistical 
errors without significant loss of structural detail and also to estab- 
lish a standard energy averaging interval (AE = 50 kev) for presenting 
measurements at many angles. 

Insofar as we were concerned with pure Wigner cusp effects, we 
could not use data averaging procedures in the experiments on lithium 
and iron. Here, tne magnitude of the structure in the cross sections 
turned out to exceed that expected ol scuspecttcers; inethe expenmnienus 
on iron, the structure was particularly marked. A successful search 


for a cusp effect would depend upon distinguishing this basic structure 
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Fig. 4. Neutron differential elastic (90°) cross section of 
natural cerium for three widths of energy averaging interval. 
(a) Original data; AE = 7 kev; statistical errors +5%. (b) Same data 
averaged over AE = 25 kev; statistical errors 2.5%. (c) Result of a 
larger averaging interval, QE = 50 kev; statistical errors approxi- 
mately 1%. The averaging in (c) was selected as the standard for cusp 
analysis of cerium. The loss of structural detail upon successive 


averaging is apparent. 
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in the cross sections from a true cusp effect. We illustrate with 
Fig. 5, which gives iron total cross section measurements, and which 
shows that use of 25 kev incident neutron energy resolution washed 
out significant structure elucidated by experiments with / kev spread 
in incident neutron energy. 

In addition to washing out structure, a poor-resolution inci- 
dent beam also has the well known effect of distorting sharp features 
in cross sections. We observed this in differential cross section 
measurements; a comparison of Fig. 11 and Fig. 14 shows a feature shift 
of about +10 kev as the result of using 30-kev compared to 7-kev energy 
resolution. This amount of feature shift would make a cusp analysis 
impossible. 

In reporting some of the lithium differential cross sections 
we have energy averaged the data for clarity in presentation (see Fig. 
10). It is clear from the foregoing discussion that in so doing we 
conceded that there was no cusp effect in the data and that remaining 
structure must be attributed strictly to the elastic scattering cross 


section of the lithium nucleus. 
Data Analysis 


Differential Cross Sections 


We give our calculation of the measured differential elastic 
cross sections using (1) the basic data, which consisted of the counts 


C = C, (scatterer in) - C. (scatterer out), (2) geometrical factors, 


0 








Fig. 5. The total neutron cross section of natural iron, using 
neutron beams of two energy spreads: (a) AB = 25 kev; (b) AE =7 
kev. The size of the symbol indicates the size of the statistical 


errors. 
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and (3) an overall efficiency measured in a calibration run. In order 
to make a detailed comparison of gross structure, cusp effects in 
cross sections, or angular distributions with theory, or with other 
experimental work, additional corrections for multiple scattering and 
finite scatterer effects had,to be applied. These corrections are dis- 
cussed in Appendix E and tabulated in Table III. 

Measured cross sections will be given in the center of mass 
(CM) frame and will be expressed as functions of mean polar CM angle 
(@ (CM) ) and mean CM incident neutron energy (E (CM)). Our reason 
for presenting CM data is that theoretical formulas for cusp phenomena 
and threshold energies are more conveniently expressed in CM terms. 
The quantities do(CM), @(CM), and E(CM) were computed from the cor- 
responding quantities expressed in the laboratory frame, using well 


oy) 


known transformation formulas. In the present experiment the differ- 
ences between lab (laboratory) and CM quantities were negligible except 
in the cases of iron and lithium. 

E(lab) was determined by subtracting one-half the incident 
neutron energy spread from the maximum lab energy of the neutron beam. 
Incident neutron energy spread was computed, assuming statistical in- 
dependence of the two contributions, from (1) the geometrical spread, 
fixed by the neutron collimation angle and T(p,n) kinematics to be 
4 kev, at 1.5-Mev neutron energy; and (2) the proton energy spread 
caused by the thickness of the metal target, containing tritium. The 


latter spread was computed from the target thickness at T(p,n) thresh- 


& 
eee [1.019-Mev proton energy (lab)] as determined from the slope of 
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the rise at threshold and conversion factors relating proton stopping 
power of the target at threshold to the stopping power at the energy 
of ‘itera Neutron energy spread is directly related to proton 
energy spread by T(p,n) ae 

Maximum incident neutron lab energy present in the incident 
beam was determined from maximum proton lab energy by the same kine- 
matical relations, where the proton energy was calibrated in two ways. 
First, by using the accurately known T(p,n) threshold energy we ef- 
fectively calibrated the Van de Graaf machine analyzing magnet at 
1.019 Mev. The magnet field was measured with a proton magnetic res- 
onance probe, between the pole faces, to obtain with high precision 
an RF frequency which could be associated with 1.019 Mev proton energy. 
We used the intersection of the straight line extrapolation of the 
lead slope of the T(p,n) yield versus frequency curve and the frequency 
axis as indicating the true threshold. Since proton resonance fre- 
quency and the magnetic field (B) are linearly related (f =p ° B/h, 
where f is frequency, UL is the proton magnetic moment, andh is 
Planck's constant), and the analyzing magnet is a momentum selector, 
then maximum proton lab energy is E, (lab) = 1.019 oe Mev where 


a 


fon is the previously described threshold proton resonance frequency. 


12 
Second, because C has a sharp resonance at 2.076-Mev lab 
neutron energy, a long counter of the type used in this experiment has 
an anomaly in its efficiency versus energy curve at 2.076 Mev due to 
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the increased cross section of C in the paraffin moderator. ious, 


a dip occurs in monitor neutron counts at this neutron energy. The 





proton resonance frequency at this dip established a calibration point 
for the mean neutron energy which served as a check on calibration by 
the first method. The two calibration points conveniently straddled 
the proton energy region involved in the present experiments. 
Initially, the "effect counts’ C were normalized as described 
in the first section of this chapter, giving normalized effect counts, 
Cy" 
For a given scattering sample and incident neutron lab energy 


BE the relation between C and the differential cross section is; 


N 
Cy =% f(z = ER) RL] f(E - E.) a 2. 
alleeecp| = nd, a)] ao(é@, E)/o,, (27) 


where Y is the neutron yield of the target per unit solid angle; 

f[(£ - En)» R is the intrinsic crystal efficiency at the lab energy 
of the scattered neutron and at the crystal-to-scatterer distance, R3 
ER is the lab recoil energy of the residual nucleus; ty 


a discriminator correction to be discussed; 2 and 2. are the solid 


(EB - E,) is 


angles subtended by the crystal to the scatterer center and by the 

scatterer to the proton beam spot, respectively; nd is the number of 
scattering nuclei per unit area in the direction of the incident beam; 
0, is the total cross section of the scatterer material; and do(@, E) 
is the differential elastic cross section per unit solid angle. Since 


do(@, E) has been assumed constant [otherwise Eq. (27) would contain 


integrals over solid angle of the finite detector and energy spread of 
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the beam] over the angular region described by Q, and the energy 
region included by the energy spread of the incident neutrons, the 
symbols @ and E refer to mean polar angle and mean neutron beam 
energy respectively [i.e., 9 (lab) and E(lab)]. 


The corresponding formula for "calibration counts’ is 


a XQ, f (E,R,) £(E), (28 ) 


where 2 is the solid angle subtended by the crystal to the proton 
beam spot and R. is the associated scatterer to tritium target dis- 
tance. 

Taking the ratio Eq. (27)/Eq. (28) yields 


f (E,R) £,(E) 0. 


do(@,E) = = TES EY RT RES 7 a a 


X Ggfl - exp(- no, 4)]"~ K(Cy/C_) (29) 


where K is a constant normalizing the effect and calibration counts 
to the same accumulated proton charge. 

The factors fy are necessary corrections which arise because 
the discriminator bias on detected linear pulses was fixed. As a re- 
sult, the counted fraction of the total detected elastic neutron spec- 
trum is a function of energy. This is shown in detail by Fig. 6 and 
is explained further in the figure caption. In the detection of scat- 
tered neutrons, the factor f£ becomes £(E = BE.) because the lab 


energy of the scattered neutrons rather than that of the incident neu- 


trons determines the discriminator correction. In all except the 
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Fig. 6. Linear pulse spectra for elastically scattered neutrons 
at three energies, showing the energy dependence of detection efficiency 
due to fixed integral discrimination bias. The spectra were obtained 
with the detector at 125° (lab) to a lithium scatterer and are the re- 
sult of subtracting the spectrum with scatterer out of the neutron beam 
from the spectrum with scatterer in the beam. This figure only illus- 
trates an experimental point in the data analysis. Actual runs were 
made with the linear bias set at a fixed voltage, corresponding to ap- 
proximately 1/2 maximum pulse height for a neutron energy equal to the 
mean energy of a complete run. In this example, the bias would be set 
at about pulse height 70, but the run could not extend as low as 0.28- 


Mev neutron energy. 
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experiment on Li and Fe f,(E - E,) = f(E) because the lab recoil 
energy of the residual nucleus was negligible. Thus the ratio 
f,(E)/f)(£ ~ E,) in Eq. (29) was ~1. For the lithium and iron experi- 
ments we determined this ratio by comparison of pulse height analyzer 
spectra of linear pulses taken during "calibration" and "effect" runs. 
We simplified the ratio f (E, R,)/£,[(E Z E,) RI, which 
is a ratio of intrinsic crystal efficiencies, in a reasonable way, by 
assuming when necessary (in lithium and iron experiments) that the 
energy dependence of Ss was separable from the distance dependence, 


R,.- Recalling the calibration equation, Eq. (28), we see that 
c /[I¥(E) 2, £,(E)] = £(E,R ) (30) 


Thus, if we divide the plot of on versus E by the incident flux, 
Y(E), obtained from monitor neutron counts, and by f,(E), obtained 
from linear pulse spectra, we see that the energy dependence of f. 

is determined. The same energy dependence, by assumption, applied to 

f (E - BoR,)s hence by translation of the curve La along the energy 


axis by an amount E_ (calculated from kinematical relations~”) we 


R 
obtained the ratio f(E - ER) /f_(E, R) insofar as the (separated) 


energy dependence is concerned, i.e., 
f(E, R)/f,((E - E,), RJ = e,(E) £(R,)/f,(R), (31) 


where e) (EB) is the just determined energy dependence of the f, ratio. 
It will be useful to write 


2/9, 2 = e,(R,“/A_)» (32) 
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where A. is the scatterer area projected perpendicularly to the beam 
axis and €) is a small correction for the degree to which the three 
solid angles are not exactly expressible as projected area/distance 
sHluared. 
Substituting (31) and (32) in Eq. (29) we obtain 
£ p=) 


R.) 
do(@,E) = e)(E) a e,(R, */h. ) FEE) 


x O,n( 1 - exp(- n On a oe K oF OE (33) 


The quantity ¢, Sa) /feCRe) = e(R) represents the inverse 
of the effective crystal efficiency as a function of distance, R 
normalized to its value at R.. = R.: If we make the reasonable assump- 
tion that this quantity is very close to unity for R,. equal to and 
greater than R | (R, was 28 in. in our experiment) then neutron counts 
versus distance can be compared to inverse square law behavior (normal- 
ized at large distances) to obtain é,(R_). This measurement was made 
using a PoBe neutron source at several distances (R.) .from the 
crystal, to obtain relative counts in unit time for each position. 
The results are shown in Fig. 7. Large corrections would be necessary 
for distances less than 5 in. In the present experiments the factor 
én{R, ) was typically 1.10 to 1.20 (R varied from 5.5 to 9 in. de- 
pending upon the polar angle @ of the cross section measurement). 

In view of the foregoing reductions we may write the most use- 


ful form of the expression for cross section: 
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Fig. 7. The deviation from We nae ae efficiency response of 
a 2" dia xX 1" long stilbene crystal, oriented with a diameter parallel 
to the distance coordinate. Curve A; Hees a fitted to the data 
with largest abscissa. Curve B: Smooth curve through the experimental 
points. Curve C;: Calculated curve using graphical analysis and effi- 
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ciency corrections from Marion and Fowler. 
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oon Gou?Boy) = (Ry /Ag) €o6R,) €y €, (EVE ,(2)/t,( - BR) 


x O,,( 2 - exp(- n Orn ain K Cy(E)/C(E) (34) 


where we have included a factor E3 denoting the conversion of the 
cross section to the CM frame. For all except the experiments on iron 
and lithium, the factor ¢) (E) f (E)/fp(B - ER) was unity, which 


simplified data analysis. 


Total Cross Sections 
The data analysis for total cross sections is described by 


0,(8(0M)) = €4(ND)"~ Inf(Cy - Cy)/(Cy ~ CQ)} (35) 


where 0,,( E(CM) ) is the total cross section as a function of mean 
center of mass neutron energy, 25 is (as before) a symbol of the 
conversion of the energy dependence from lab to CM frame, ND is the 
number of nuclei per unit area in the sample, ln means natural log- 
arithm, Ch is the number of direct neutron counts, Ce is the number 
of background counts, Cr is the number of transmitted counts, and 
where Cy Cy» and Ch are all functions of neutron energy and are 
all normalized to a fixed accumulated proton charge. No monitor neu- 


tron normalization was applied because only count ratios are involved 


with the machine energy at one setting. 


~ 
- 





at 


Error Analysis 


Differential Elastic Cross Sections 


Contributions to the uncertainty in the absolute magnitude of 
the differential elastic cross sections were present in each factor of 


Eq. (34), reproduced below: 


do, v6 Pou? = n(cm)? . (R“/A,) e,(R, ) Ey e) (E) (f,(E)/t)(£ = E,)) 


x 0,.( 1 “Jers ie K Cee (34) 


Statistical errors entered the results through the factor 
C,/C,- Those in on were negianibile. Tose whichemere in Cw are 
separately described in the figure captions for each set of experimental 
data. In general, the statistical errors were rather small (from 1 to 
s percent); hence, they had little effect on the precision of the abso- 
lute differential cross section measurements. 

Experimental uncertainties in the factors of Eq. (34) are 
discussed below: 

La geac Gor Cy/C There were two experimental uncertainties 
in C,. The detector was located at the scatterer "in beam" position 
during the calibration run. Its distance from the neutron source was 
about 28 in. The precision with which the crystal could be placed at 
the exact scatterer position was nat ite: in. Hence, there was an +1 
percent uncertainty from this source. Also, we include in C. an. 


error associated with the precision to which the crystal was oriented 
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with respect to the incident neutrons. The aspect of the crystal rela- 
tive to the scattered neutrons, when the crystal was in detection 
position (see Fig. 1), should have been the same, because deviations 
produce differences in effective crystal efficiency. We assign an ad- 


ditional uncertainty of ft2 percent to C as a measure of the preci- 


N? 
Sion with which the crystal aspects were adjusted. 

2. Factor (R“/a,). The experimental uncertainty in R, 
the crystal to scatterer distance, was +1/8 in. for mean R, about 7 
insee thus, R.. had an extreme uncertainty of 12.5 percent. The dimen- 
sions and aspect (with respect to the beam axis) of the scatterers were 
not precisely known, which led to a 11 percent uncertainty in A: A 
total +6 percent uncertainty was assigned factor R./A. 


-1 
3. aebaAeitor 6,,( 1 - exp(- n On d)) . Uncertainties in the value 


of o enter this factor only in second order terms in no 


- d. The 


T 
precision to which the factor nd was known determined the experimental 
error. A tl percent uncertainty was assigned this factor. 

4, Factor K. Although the current integrator used in deter- 
mining the normalization to fixed proton charge is rated at +1 percent 
reproducibility, the mean normalization for several counting periods 
of the same accumulated charge was quite precise. So also was the mean 
normalization for several counting periods of a different fixed proton 
charge (i.e., the counting period for calibration points was much less 


than that for runs involving detection of scattered neutrons). The 


uncertainty in K was at most t0.5 percent. 
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5. Factor e,(R). This factor was obtained by comparing ex- 
perimental points with a theoretical (atewnmmn) behavior for effective 
efficiency versus distance (see Fig. 7). Figure 7 shows that the data 
did not precisely fix the normalization of the curve at large distances. 
The figure also shows some discrepancy between experimental and calcu- 
lated points (discussed elsewhere in this report). Considering = 7 
in. typical, we estimated +3 percent uncertainty in e,(R.). 

6. Factor E.: This factor is symbolic, representing conver- 
sions involved in quoting results in terms of mean CM energies, CM 
angles, and CM intensities, where CM is center of mass. Quantities 
M, /M, and M, /(M, “ M,), which are precise in the present usage, 
dominate these conversion formulas, the uncertainty in lab angles and 
energies playing a minor role. An extreme uncertainty +0.5 percent 
was assigned factor Sas 

7. Factor e) (E) f,(E)/f£)(B - BE.) In all except the lithium 
and iron experiments, this factor is unity to good precision. In 
lithium and iron cross sections it is a source of additional uncertainty. 
These uncertainties, which are potentially large for a single data 
point, are reduced if one observes a smooth trend in the factor when it 
is determined for points at many energies. We estimated a +2 percent 
uncertainty for iron cross sections and +5 percent uncertainty for 
lithium cross sections from this factor. 

some of the experimental errors, given in 1 through 7, were 


2 
correlated; for example, R, is anti-correlated to e,(R,). For the 


purposes of this experiment, detailed examination of such correlations 
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was not necessary and the procedure adopted was to assume independence 

of the errors; thus, the sum of the above percent uncertainties applied. 
We concluded that exclusive of statistical errors and those 

errors associated with multiple scattering corrections (see Appendix 

E), the absolute differential cross sections of zirconium, barium, and 

cerium were known to within +14 percent; of iron to within +16 percent; 


and of lithium to within +19 percent. 


Total Cross Sections 


Measurements of total cross sections were accomplished in a 
neutron background of less than 5 percent of direct neutron counts which 
minimized second order background corrections to the data. The samples 
had small dimensions perpendicular to the incident beam direction, be- 
ing just sufficient to completely shield the detector from direct 
neutrons. The sample-to-detector and sample-to-source distances were 
large compared to the detection crystal dimensions. Neutron transmis- 
sion through the samples was approximately 60 percent of the incident 
flux. All these are excellent experimental conditions for measurement 
of total cross peerionk A crude calculation indicated that in- 
scattered neutron corrections to the cross section were at most 2 per- 
cent. The main constraint on the amount of in-scattering was the 
small solid angle the detector crystal presented to in-scattered neu- 
trons. An uncertainty of +3 percent was assigned to account for second 
order background scattering effects and sample in-scattering corrections 


which we did not calculate in detail. Thus, we are assuming in effect 
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that the uncertainty in our rough calculation is larger than the cor- 
rection. A +l percent extreme error was estimated for the measurement 
of ND, the number of scattering nuclei per unit area. 

Hence, exclusive of statistical errors which are given sepa- 
rately for each experiment, the total uncertainty in our absolute total 


cross section measurements was +4 percent. 
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CHAPTER V 
EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 


Energy and Angle Dependence of the Cross Sections 


Although the basic purpose of the experiment was to search for 
threshold effects in neutron cross sections, data such as angular dis- 
tributions and differential elastic and total cross sections were 
obtained which are interesting in their own right, exclusive of possible 
cusp effects. These data will be discussed briefly, with emphasis on 
illustrations to clarify the discussions. 

The measured angular distribution of 1.6-Mev neutrons scattered 
elastically from natural cerium is presented in Fig. 8. The data 
points are cross plotted from measurements as a function of energy at 
separate angles. As shown in the figure, an optical model prediction 
Gis Moldauer,~' modified for multiple scattering and finite scatterer 
effects (see Appendix E) agrees very well in shape and, except at back 
angles, in magnitude, with this data. 

Figure 9 shows the measured angular distribution of 600 kev 
(CM) neutrons elastically scattered from natural lithium. The data of 
this experiment compares favorably with that from a compilation of 
other experiments by Hone eonat which is also shown. ‘This substan- 
tiates our estimate that multiple scattering corrections were small for 


the lithium experiments (see Appendix E). The shape of the distribution 
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Fig. 8. Angular distribution of neutrons elastically scattered 


ak 
t optical 


from natural cerium at E| = 1.6 Mev. Curve 1: A Moldauer 
model prediction. Curve 2; Curve 1 modified with a multiple scattering 
correction. Curve 3: Curve 2 modified with calculated experimental 


finite scatterer effects. Triangles are experimental points from the 


present experiment. 
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Fig. 9. Angular distribution of neutrons elastically scattered 
from natural lithium at E = 0.6 Mev. Open circles are data from the 
present experiment. This data is not corrected for multiple scattering. 


= 
Triangles are data from Howerton's compilation. 
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is interesting, being peaked backwards, presumably as a result of inter- 
ference between the strong p wave nasa ee at 0.25 Mev and an s 

wave background. This backward peaking means that one should look at 
forward angles for cusp effects, because the theoretical percentage ef- 
fect is inversely proportional to (aa)1/2 (see Chapter II). 

In Fig. 10 we show some of the differential cross section data, 
versus energy, from which the angular distribution, Fig. 9, was cross- 
plotted. At the lower energies the upper edge of the p wave resonance 
at 0.25 Mev is in evidence. The wavy structure in all the cross sec- 
tions is probably not an experimental effect; it is quite surprising, 
because for a light nucleus such as lithium one might expect rather 
smoothly varying cross sections at these energies, except for isolated 
resonances. The inelastic cross section to the first excited state of 
te as measured by Freeman, Lane, and Rese? shows similar structure. 

The differential elastic cross section of iron, versus energy, 
was measured for four angles using a rather thick neutron target (#3, 
Table I) in the vicinity of 0.847 Mev. The results are shown in Fig. 
11. It was later determined that a much thinner neutron target was 
needed to resolve the structure and to search for cusps (see Fig. 14; 


target #5, Table I was used). 


Energy and Angle Dependence of Threshold Features 


We discuss here the experimental results as they pertain to 
Wigner and "energy-averaged" Wigner cusps. Frequent reference will be 
made, explicitly and implicitly, to the predictions of cusp theory, 


Chapter II. 
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Fig. 10. Neutron differential elastic cross sections (CM) of 
natural lithium versus center of mass (CM) energy for five angles (CM). 
The data are not corrected for multiple scattering. The original 
experimental results have been averaged such that statistical errors 


in these points are approximately +1%. 
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Fig. ll. Neutron differential elastic cross sections of natural 
iron versus energy for four angles. The data is not corrected for 
multiple scattering. Incident neutron energy spread is approximately 


25 kev. Statistical errors are indicated by the size of the symbol. 
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Lithium (92 .6% ri!) 


The cusp results for lithium were essentially negative. The 
theory of Appendix B, Case I, applies. We saw in the preceding section 
that the largest percentage cusp effect in dag is expected at forward 
angles. Figure 12(b) shows that an experimental effect was indeed 
present at 34° (CM) in the vicinity of 0.478 Mev, the first excited 


t 


state of Li. However, using the reaction cross section magnitude of 


e 
7 and the applicable theory, we find that the 


Freeman, Lane, and Rose, 
observed effect was at least a factor 5 too large to be explained as 

a Wigner cusp; an effect of approximately 2 percent is predicted 50 kev 
from threshold. Figure 12(a) compares our measurement of lithium total 
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cross section with that from Howerton's compilation. Although the 
two sets of data agree, no Wigner cusp is in evidence in either. We 
had hoped to see an effect in the total cross section, because the re- 
sults of Lane, et al., (see Fig. 13) appear to show a break near 0.5 
Mev in the coefficient BY of the cross section expansion (see Chapter 


II). The energy scale change in Fig. 13 at about 0.5 Mev may bias this 


interpretation of their data. 


Iron (91.7% Fe’) 


With iron distinct features appear in all cross sections at the 


6 0.847 Mev). This is 


+ 
energy of interest (2 first excited state Fe 
exemplified in Fig. 14 where a resonance feature appears, apparently 


at an identical energy (just above 0.847 Mev) in the total cross sec- 


tion, in the 145° differential elastic cross section and in the reaction 
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Fig. l2. Neutron cross sections of natural lithium versus 
energy. (a) Total neutron cross section; triangles are data from 
Howerton 's aemeiAtion: © and solid circles are the energy-averaged 
(A E = 20 kev) data from the present experiment. (b) 34° CM differ- 


ential elastic cross section of the present experiment. 
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Fig. 13. Differential cross section for scattering of neutrons 


by Lil. Reproduced from the work of A. Langsdorf, Jr., R. O. Lane, and 


2 
J. E. Monahan. o 
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Fig. 14. (a) Neutron differential elastic (145°) cross section 
of natural iron versus energy. (b) Total neutron cross section of 
natural iron. (c) Reaction cross section Feo” (n, n', y (.847 mtev)) Fe”, 
Statistical errors are the size of the symbols in (a) and (b), but are 
+10% in (c). The smooth curves show the result of removing a calculated 


Wigner cusp effect under one possible assumption on the phase of the 


scattering amplitude. 
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eross section ea, n')Fe® (0.847 Mev). The latter cross section 


ee 


was measured by Tucker, et al., using the same tritium target. 
other resonance feature which seems to be a ‘mirror image’ of the 
aforementioned resonance across threshold is evident in these data, 

and in all the differential cross sections at other angles investi- 
gated [see Fig. 15 (six angles)]. Application of the most general 
theory is hopeless because of the large number of unknown phases. At 
least six phases contribute to the cusp, even with a restriction on 
participating incident partial waves to those with orbital angular 
momentum £ less than 3; see Appendix B. A reasonable and tractable 
special case of the theory, which might apply to the iron data is dis- 
cussed and explicitly evaluated in Appendix B, Case II. Unfortunately, 
the special assumptions involved allow cusp analysis only over a nar- 
row energy region, entirely above threshold. The results of this 
analysis are superimposed on the data in Fig. 15 for differential cross 
sections (6 angles) and in Fig. 16 for the total cross section. The 
results are neither startling nor conclusive, although by removal of 
the cusp effect the peak energies and widths of the above threshold 
resonance in the differential elastic, total and inelastic cross sec- 


tions can be brought into exact coincidence. 


Cerium-Barium-Zirconium 


We consider now the results obtained for targets of higher 
atomic number in the region where the optical model is expected to ap- 


ply because compound resonances are closely spaced--cerium, barium, 





Fig. 15. Neutron differential elastic cross sections of 
natural iron versus energy for six angles. The data is not corrected 
for multiple scattering. (a) gives the data over the full energy 
range of the run. The lined box shows the section of the data which 
is expanded into (b). In (b) statistical errors are indicated by 
the symbol size. Open triangles show the result of removing a cal- 
culated Wigner cusp effect, under one possible assumption on the 


phase of the scattering amplitude. 
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Fig. 16. Total cross section of natural iron. The size of 
statistical errors is indicated by the size of the symbol. Triangles 
show the result of removing a calculated Wigner cusp effect under one 


possible assumption on the phase of the scattering amplitude. 
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and zirconium. The theory of "energy-averaged " Wigner cusps is then 
also expected to apply (see Chapter II). 

Composite data for cerium (88.5% get”), given in Fig. 17, is 
representative of the cusp analyses we present for all three elements. 
A special feature of this figure, however, is that we present the 
original data, to show the degree of scatter on interleaving runs (most 
of which is not statistical, but experimentally reproducible). In the 
remaining figures presented for this group, the data was averaged over 
an energy interval suitable to reduce statistical errors to about tl 
percent (refer to Chapter IV). This procedure standardized data presen- 
tation for a large number of angles. 

This and other figures in this group also show cross section 
trends predicted by the optical model. These were obtained from the 
work of Campbell, et Pad The optical model shapes were interpolated 
from very few points given in the reference and should not be consid- 
ered a precise computation; hence, they were normalized to the experi- 
mental data at threshold energy. The magnitudes of the predicted cross 
poomiens-. were consistent with the estimated experimental errors 
(Chapter IV) and the corrections given in Appendix E. 

A cusp effect is quite apparent in Fig. 17(b) by comparison of 
the optical model trends with the experimental cross section. To com- 
pare the magnitude of the cusp effect with theory one could subtract 
the experimental cross section from the optical model prediction and 
plot the difference as a function of neutron energy. But, as mentioned 


before, the optical model trends are somewhat uncertain and, furthermore, 








Fig. 17. The experimental inputs to a cerium "energy-averaged 
Wigner cusp analysis. (a) Energy-averaged (AE = 40 kev) experimental 


inelastic cross section cet O(n, n' ¥ (1.597 Mev) ) gett? 19 (bd) 


Origi- 
nal data for elastic 95° scattering of neutrons from natural cerium. 

The two symbols indicate the reproducibility of the data on successive 
interleaving energy runs. The size of the symbols indicates the size 
of the statistical errors. (c) Same data as (b) corrected, point by 
point, for a calculated "energy-averaged”" Wigner cusp. The smooth 
curves in (b) and (c) are optical model predictions of Campbell, et al., 
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normalized at 1.6 Mev. 
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the experimental cross section has natural fluctuations, which will be 
discussed below. Hence, it was thought more appropriate to remove from 
the experimental cross sections the theoretical cusp effect and to show 
(a) that in the resultant curve there would be no more cusp effect--at 
least within the natural fluctuations of the cross section, and (b) 
that the resultant curve would follow the optical model trend, which 
theoretically represents the cross section without threshold effects. 
The procedure for removing a theoretical cusp from the measured 
cross section was to multiply the reaction cross section (i.e., inelas- 
tic neutron scattering cross section) by the isotopic percentage of the 
isotope of interest, by the appropriate multiple scattering and finite 
scatterer correction (the reciprocal of Table III (Appendix E) values), 
and by the theoretical angle dependent factor (Chapter II). This neg- 
ative quantity was subtracted from the experimental data [see, for 


example, Eq. (26)]. 


1h0 


Cerium (88.5% Ce — ) 


We now Consider the specific case of cerium, for which Eq. (26) 
applies. Tucker, et al., measured the reaction cross section [Fig. 
iy a,).) Cae (oe n')cer*( 1.60 Renae As previously mentioned, inci- 
dent d waves dominate the reaction to at least 150 kev above threshold, 
hence £' =2 in Eq. (26). The results of removing the theoretical 
cusp effect are shown in Fig. 17(c). ‘There is little doubt that the 
theoretical "energy-averaged' Wigner cusp, with £' =2 is a satis- 


factory description of the change at threshold. Note Eq. (26) predicts 





S16. 


that if incident s waves were imagined to dominate this reaction, rather 
than the actual d waves, the resulting expected cusp, at polar angle 
90°, would be a factor 20/13 or about 50 percent larger. The plot cor- 
responding to Fig. 17(c), corrected for an expected (£' = 0) cusp, 
would show an upward break at threshold and would follow neither the 
below threshold data trend nor the optical model trend. Thus we can 
conclude the reaction is not dominated by incident s waves. This type 
of information can be useful in assigning excited state spins, for 
example. 

The differential elastic cross section of cerium was measured 
at ten angles, versus energy. These results are shown in Fig. 18. 
The! points obtained after removing a calculated cusp are also given 
for each angle. In all cases, the adjusted data so obtained follow 
the trend of the below threshold data and the fluctuations in these 
is the same above and below threshold. In general, the optical model 
curves agree nicely with the modified cross sections; that is, they 
follow what the theory defines as “cross section without reaction." 
There is no reason to expect this agreement to be perfect, since the 
optical model was constructed to explain gross structures only, over 
a wide range of atomic number and energy. We should recall the assump- 


tion, made for Eq. (26), that the fraction (called R in Chapter IT) 


15 
of the reaction probability for proceeding from compound nucleus states 
of a particular J' is the same as the fraction of the dominant inci- 
dent partial wave states (all with the same orbital angular momentum 


£') which have J =J'. Other assumptions for R lead to cusp 


le 








9k 


Fig. 18. The neutron differential elastic cross section of 
natural cerium versus energy for ten angles. The data is not corrected 
for multiple scattering. Statistical errors are 41%. The points are 
an energy average (AE = 50 kev) of the original experimental data. 

For each angle the left side is the measured cross section. The right 
side is the cross section with a calculated "energy-averaged Wigner 
cusp removed. The calculation is based on the excitation of the 1.6 
Mev state in cer. No cusp contributions from the o* (1.90 Mev) or 
7 (2.08 Mev) states are included. The smooth curves are optical . 


model predictions of Campbell, et al., normalized at 1.6 fe 
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magnitude predictions differing from those used here by as much as 30 
percent, at certain angles. Although not the case in the present work, 
it is possible that a more refined experiment could be used to measure 
Rye such a measurement could assist in the understanding of neutron 
nucleus interaction mechanisms, particularly if it turned out that a 
particular reaction preferred one value of compound nucleus Je 


another, even though resonances were closely spaced and possibly over- 


lapping. 


Barium (71.7% pal38) 


The measured cross sections of barium and the results of the 
appropriate cusp analysis are compiled in Fig. 19. The reaction cross 


1380 Ainge" 19 


section Ba (1.427 Mev) was measured by Tucker, et al. 
Again the spin of the first excited 1.427 Mev level is 2” and the ground 
state 0”. That s wave neutrons dominate the reaction very far from 
threshold is not so clear in this case, since the reaction shape is 
distorted by resonance structure. However, optical model calculations 
show that barium is near a maximum in the s wave strength Functions’” 

As in cerium then, the reaction is considered to proceed by the 2£' =2 
incident partial wave. The results of theoretical comparison are sat- 
isfactory. The cusp effect is smaller than in cerium because of the 
smaller isotopic percentage and the smaller reaction cross section near 
threshold. In Fig. 19(c) we see that it would be difficult to distin- 


guish a cusp at all were it not for a comparison with the optical model 


trend which happens to be rising. 
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Fiz. 19. Weutron differential elastic (90° and 125°) cross 
‘ection. of natural bariom versus energy. (a) Energy-averaged (AE = 
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Zirconium (51.5% eee 


The measured cross sections of zirconium and the results of 
cusp analysis are compiled in Fig. 20. The inelastic reaction 
a” (ip a! eee lee Mev) was! measured by Tucker, et al.” The spin 
of the 1.75-Mev first excited state is o* and the ground state 0”, 
Tucker also showed that p wave outgoing neutrons dominate the! reaction 
not far above threshold. Hence we assume that the reaction proceeds 


90 


from incident p waves, i.e., #' =1. For Zr” the low isotopic per- 
centage and small reaction cross section near threshold combine to 
produce a small calculated cusp effect. At 60° polar angle the pre- 
dicted cusp effect is too small to be distinguished in the data. Note 
also that at 60° the optical model trend does not agree well with the 


measured cross section slope. At 130° the cusp is distinguishable and 


the theoretical cusp fit 1s satisfactory. 


Total Cross Sections 


The theoretical "energy-averaged " Wigner cusps which were used 
to explain the preceding experimental results (Ce, Ba, Zr) are based 
on the same assumptions that underlie the optical Wes Sao These as- 
sumptions also lead to predictions that the energy average of the total 
cross sections should be smooth across thresholds [see Eq. (23)]. Our 
measurements of the total cross sections of cerium and barium are given 
in Fig. 2l. The predicted cusp effect for the total cross section of 
zirconium--even if the energy-averaged theory would not apply--would 


be so small that we did not investigate this case. For completeness, 
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Fig. 20. Neutron differential elastic (60° and 130°) cross 
sections of natural zirconium versus energy. (a) is a 40 kev average 
of the reaction cross section Wen a) ae (lege Nees Statis- 
tical errors are +5%. The left sides of (b) and (c) are an energy 
average (AE = 50 kev) of the measured cross sections. This data is 
not corrected for multiple scattering. Statistical errors are +1%, 
The right sides of (b) and (c) show the results of removing calculated 
"energy-averaged Wigner cusps. The smooth curves are optical model 


30 


predictions of Campbell, et al., normalized at 1.75 Mev. 
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Fig. 21. Neutron total cross sections for natural barium, (a); 
ceriun, (b); and zirconium, (c). Statistical errors are indicated by 
the size of the symbols. The data on barium and cerium were measured 
in the present experiment. The data on zirconium is taken from 


27 


Howerton's compilation. In (b) circles are the total cross section 
of cerium; triangles are integrated elastic cross sections, corrected 
for multiple scattering; and squares are the sum of the integrated 


elastic and of the reaction cross section (corrected for isotopic per- 


centage) Ge. ln, eae gl allel ee) Noes 
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the zirconium total cross section data from Howerton is also presented 
in Fig. 21, showing the inadequacy of widely spaced energy points for 
the present study. The data for barium and cerium were taken with an 
energy spread of about 30 kev, yet they are certainly not smooth func- 
tions of energy, displaying marked structure. The mean values of these 
fluctuations do not seem to drop markedly as threshold is crossed, sug- 
gesting agreement with the theory. If, however, the theory of energy- 
averaged cusps does not apply, but one assumes instead a pure Wigner 
cusp phenomenon, the cusp predicted is at most a 5 percent effect which 
would be masked by the fluctuations. Hence, the total cross section 
data cannot be considered a verification of either assumption. 

That the fluctuations are not experimental is supported by the 
differential cross section data for cerium at forward angles (Fig. 18) 
which shows the same structure. Note that the structure does not con- 
tain recognizable resonance shapes. It is possible that the consider- 
ations of Ericson on statistical fluctuations in cross sections might 


ol 


apply here. He shows that if a compound system is characterized by 
many overlapping resonances, each of width IT, and with relative ran- 
dom phase of the value quantities TY), of resonance eee ee cross 
sections would have a fluctuating structure of average width, I. 

In the case of cerium only, sufficient differential elastic 
data was taken to permit integration over all solid angle. The total 
elastic cross sections so obtained for several energies, after correc- 


tion for multiple scattering, are plotted in Fig. 21(b). ‘The drop at 


threshold (1.6 Mev) is evident, as, in accord with any theory of 
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threshold effects, flux conservation alone requires the total elastic 
cross section to drop above threshold for a reaction. In the case of 
see the experimental differential cross sections all drop, so the 
total elastic cross section does also. It is reassuring that, when 
the inelastic scattering cross section is corrected for isotopic per- 
centage and added to the integrated elastic cross section, the result- 
ant cross section (i.e., the total cross section) is smooth across 
threshold. The agreement in magnitude between the latter curve and 
the directly measured total cross section is considered somewhat 


fortul touse 
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CHAPTER VI 
CONC LUS IONS 


The experiments have shown that under suitable conditions the 
threshold effects expected in neutron elastic scattering can be ob- 
served. From an experimental point of view particular attention must 
be paid to the choice of energy resolution consistent with observing 
a cusp effect in reasonable time. 

Explicit theoretical expressions have been obtained for ex- 
pected Wigner and "energy-averaged” Wigner cusp effects under general 
angular momentum conditions, within some reasonable restrictions imposed 
in the derivation. The effects were found to be contained in the coef- 


ficients, B of the Blatt-Biedenharn formulation of the general 


L’ 
scattering problem. Since many survey experiments involving neutron 
scattering have been analyzed in terms of the BL coeiiierents, we 
make the obvious suggestion that such results be examined near excited 
state energies for possible tests of cusp theory. 

The expected Wigner cusp in lithium cross sections was expected 
to be an approximately e percent effect, 50 kev from threshold. Our 
experiment was not sufficiently refined to pick out an effect of this 
size. 

Resonance structure in the measured cross sections of iron made 


a Wigner cusp interpretation ambiguous. Features occurred at the thresh- 


old energy, but their analysis in terms of cusp theory was inconclusive. 
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The measured differential elastic cross sections of Zr, Ba, 
and Ce dropped at inelastic thresholds. Analysis of the drops in 
terms of "energy-averaged | Wigner cusp theory resulted in good agree- 
ment between theory and experiment. This indicates that the usual 
energy-averaging assumptions underlying the optical model are indeed 
applicable to these nuclei. Also, the cusp theory ease that there 
cannot be a large amount of direct interaction between incident neu- 
trons and these nuclei at the neutron energies used in these experi- 


ments (S 2 Mev), or else below threshold effects would also be seen. 





APPENDIX A 


THE PHYSICAL BASIS FOR THE S-MATRIX AND ITS USE 


IN CROSS SECTION FORMULAS 


We will review the principles of scattering analysis and the 
definition of the S-matrix as presented by Blatt and tL ogigeiersa = 
They pointed out that general expressions for cross sections in terms 
of the S-matrix are rarely given in standard texts. 

The general formulas are applicable to any collision process 
in which two particles collide and two particles emerge. However, in 
this appendix, and after a certain stage in the development, only the 
special case of elastic scattering with no £ change will be con- 
Smeered . 

The channel index @ defines the type of incoming particle and 
the state of the struck nucleus. The channel spin, s, is the total 
Spin angular momentum of the channel, obtained by vector addition of 
the intrinsic spin, i, of the incident particle and the spin, I, 
of the struck nucleus. The orbital angular momentum (in the center of 
mass system), £, is combined with s to form the total angular 
momentum, J. Only in this particular appendix, primes will be used 
to refer to outgoing or final state quantities. J is conserved during 
the collision, hence the values of £' and s' must be consistent 
with J' =J. Similarly, the parity xn, of the initial system, is 


conserved, and a =1' 
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The probability amplitude sits Siaeees for a collision with 
total angular momentum J, from channel Qs to channel a's'£' is 
referred to as an element of the scattering matrix. The scattering 
matrix, or S-matrix, is defined by considering asymptotic forms of the 
collision wave functions, as follows. 

Consider one definite value of J and je = M. At any particu- 
lar energy, E, the channel wave number, ky and relative speed, 
Vy? in each channel Q are given by the energetics of the reaction. 
Only those channels which are "open’ (Ky real) are considered. Let 
94... be the product wave function of the target nucleus and incident 
particle. The spin and angle dependence of a wave function of total 


angular momentum J, J, = M, orbital angular momentum £, and chan- 


nel spin s is given by 


S 


= > (ésm, m_|2s Oy e.o) x (Al) 
PT. mn, ae m.=-s £ | £m , Sm. 


where *s ym, is the incident spin wave function and Lon is the well 
known spherical harmonic of angular momentum £ and projection nm. 
Using these definitions, the most general wave function in channel @Q,s 
with total angular momentum quantum numbers J, M, consists of the 
superposition of incoming and outgoing spherical waves. Asymptotically 
(at sufficiently large i stan we can write (for neutral particles) 


Vo men ere Qs 


JM ; 
- Boog expl +ileyry 


JM : i: 
x {fase expl-i(k r, = 5 é x) | 


: L x1} (A2) 





ae 


Substitution of Eq. (A2) into the wave equation yields a unique connec- 


tion between the outgoing wave coefficients, B, with the ingoing wave 
coefficients, A. Thus 


JM af JM 
Bay's 2! 7 = 2 A Pa's'£' sas 0 Aas £ (A3) 


Eq. (A3) defines the scattering matrix, §. With these definitions it 


is evident that S is unitary (conserving the wave function normaliza- 


tion in time), symmetric (reciprocity theorem), and independent of M. 


It remains to decompose the incident plane wave into the func- 


tional forms we are using. The incident wave is 9, exp(ikz x, 


First, one writes 


Ss 
ae 
exp(ikz)X, =(ue)e 2 i°(22 +1) é J oy p (kr) Yh Os on 
Ss £=0 S 
where 
L+s J 
Te a = > by 


M 
(£som_| £s oH) 
s  J=|é-s| M=-J i mee 


Using the asymptotic form of Jo, 2 for kr>> ££ the following ex- 
2 
pression is obtained for the plane wave in channel a,s 


with spin 
direction given by m_:; 





1 

: s ‘Jj J+s 
a eles) = = 
exp( ik 2q)%q Pas We Pas 2 2 2 


=e J=-0 M=-J £=|J-s| 


x 


_f 5 M 
(4som_| fs JM) i°(2£ + 1) ane (AL) 


x fexpl -i(k yr, -32n)] - exp[ti(k r - 3 2 x) 
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Comparing Eq. (A4) with the standard form, Eq. (A2), we see that the 


amplitudes of the ingoing spherical waves are given by 


JM ae L i 
A = ix _(x ey) (£som, | £s JM) i°(24 + 1) 


asl a (A5) 


The amplitudes of the outgoing spherical waves are determined from Kq. 
(A3). As in the main text, we now specialize to elastic scattering 


with no change of £ or s. They are 


J+S 
JM a 
B d Fee. Ve OE 
As £ a a t=|J-s| 
x (£som_|£s JM) i*(2¢ + 1)? S 4 (A6 ) 
S ask;ask 


In order to investigate cross sections we only need that part of the 
plane wave solution, Eq. (A4), which is due to the scattering interac- 


tion; 1.e., we want to obtain y from y =y 


scattered incident 


y Therefore, we substitute Eqs. (A5) and (A6) into Eq. (A2) 


scattered 


and subtract the incident wave, Eq. (A4); the result is 


a 1 co Al J+s 
scattered a ty J Oued £=|J-s| 
ele 
2 


x (fsom., | £s JM) i*(22 + 1) 


, nl J M 
. explilkyry 2 4x)] x [1 - See nseat hits (AT) 


: MC, : 
If the spin angle function "Y., | is decomposed using Eq. (Al), then 


Eq. (A7) can be used to define the scattered amplitude, ya 
s? S 





_ exp ikr 
¥ gq osm )=i ae 2 Qs 
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where q is explicitly 


co J Jt+s £ 
—— & y 
J=O0 M=-J £=|J-s| p=-2 


Iysm'! asm (8,9) 
5? Ss 
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m(2L + 1)*( £som_| 2s JM) ( £sum’ | £s JM) 
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The differential cross section for the scattering sm. ~ dsm. is 


given by 


a (8,9)|° an (A10) 


2 
"osm'sasm_ ~ Koy bisa ;asm 
S S S 
For unpolarized incident beams, we must average Eq. (Al0) over initial 
spin direction m_ and sum over final spin direction ma: 


(A11) 


|. 
= a3 
do ; (2s + 1) 2 


Ss 
2 
=e 


Oo 
asm‘ ;dsm 
s s s 


=S ul 
Ss Ss 


This expression is greatly simplified by sum rules for Clebsch-Gordan 
coefficients, and it is from this calculation, as performed by Blatt 
and Biedenharn, that one obtains the useful final cross section expres- 
sions, involving Z coefficients. 

Finally, to obtain the differential elastic cross section, one 


must average over the possible values of incident channel spin, s. 


en 2s +1 


Boy = saliaaad (21 + 1)(2i + 1) *"as;as (Al2) 


Equations (Al2) and (All), as explicitly worked out in Ref. 12, are 


just Eqs. (ljge(2) and (3) of the mann text. 


APPENDIX B 
SPECIAL CASES OF INTEREST TO THE PRESENT EXPERIMENT 


Case I. Cusp effects in neutron scattering from Li 
near the threshold for excitation of the first 
excited (0.478 Mev) 1/2 state of tl eon 
state spin and parity, 3/2 ). 
The large resonance in the elastic cross section at about 0.25- 
+ 
Mev neutron energy has been analyzed and found to correspond to a 3 
+ 
state in 148 28 The reaction cross section a Lay ai yea | (0.478 Mev) 
has a broad resonance at 1.35-Mev neutron energy, probably associated 


eg 


with a compound 1 state. Also, the inelastic scattering is domi- 
nated near threshold by outgoing s waves, with compound angular momen- 
tum a = ie Finally, any process involving incident partial waves 
with £>1 is unlikely at 0.5 Mev because kKR=0.3. We conclude 
that in the vicinity of 0.48-Mev incident (center of mass) neutron 
energy only s and p waves participate significantly in elastic scatter- 
ing. Also, incident s waves dominate the inelastic scattering. The 
latter conclusion results from the known spins and parities of the 
ground and first excited state of ma and the assumption that no 
process with £>1 is likely. With s waves dominant in the outgoing 
reaction channel the only possible incident channel leading to the re- 


action is characterized by J' =1, s'’ =1, and #£' =0. Using these 


values, and I = 3/25 ie /2, in Eq. (12), one obtains 
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Restricting ourselves to £< 2, as indicated above, we need these 


Z coefficients 


2 
Z°(0101, 10) = 3, z°(1101,11) =o, 2204 PONT) = 5 
and all other 2° ( £501, 11) values are zero. Therefore, 
a oe Se 
el OL 
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1 1 di 
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Equations (B2) might be compared to data such as that of Langsdorf, 
Lane and isaanen, reproduced herein as Fig. 13. Unfortunately, the 
energy scale in this figure changes in the region of interest (0.5 Mev 
lab). We present the figure primarily to illustrate the utility of 
cusp expressions in the form of Eqs. (B2). The cusp prediction for 
the differential elastic cross section is obtained by substituting Eq. 


(Be) in Eqs. (9-11) to obtain 
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Using Eq. (16), we obtain the total cross section cusp term, 


dk 


cos 25 
Ao, = (3/8) (poenp 4 - } (Bh) 


; EB 
- sin e551 


where in Eqs. (B3) and (BL ) is the experimental reaction 


(91.90 exp 
cross section near threshold. 

For the casesei taf we recall the restriction imposed on the 
derivation leading to Eq. (12), that the amount of scattering with 
change of s, but no change in energy or orbital angular momentum {f, 
be small. There are no good reasons for believing this to be true in 
the present case. Change of 3 =I+i from js| =2 to Is, = 1 
(Ss. means outgoing ‘elastic" channel spin), or vice versa, in the ri! 
ease is essentially a spin flip folllowed by a recoupling of s and 
£221), and this is just the type of interaction involved in the 
reaction leading to the first excited state. We have seen, however, 
that in fact s waves (£' = 0) dominate the ti! er, thus "spin flip’ 


elastic scattering (which requires £2 1) does not couple to the reac- 


tion. Hence, Eqs. (B3) and (B4) apply. Using the reaction cross section 





La 


2 
7 and the 


value (100 kev above threshold) of Freeman, Lane, and Rose, 
present experimental values of do (3) CM), we can estimate the ex- 


pected percent effect to be 


A do/do® = (3/32 x) o_,/ao° =~ (3/32 x) .05/.05 = 3 percent effect 
(~ 100 kev from threshold) 


Correspondingly, for the expected cusp in On we obtain 


IN gic = (3/8)(.05/1.1) = 2 percent effect 


(100 kev from threshold) 


Case II. Cusp effects in neutron scattering from 
iron near the threshold 10r excrcario0n, Of sene 
first excited state (2°, 0.847 Mev) of Fe? 
(ground state spin and parity 0°). 

Equations (9) to (12) for pure Wigner cusp effects apply to 
this case because the resonance structure in the cross sections is ap- 
parently resolved. We must compare Eq. (9), the expression for change 
in cross section, with experiment. There are too many possible phases 
involved in a general treatment of Eq. (9) to make cusp interpretation 
meaningful; therefore, we sought a tractable and reasonable simplifica- 
clon . 

The experimental data, see Fig. 14, display a resonance at an 
identical energy in the reaction cross section, in the total cross sec= 
tion, and at six angles in the differential elastic cross section. On 


this evidence, we assumed that in the vicinity of this resonance the 
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scattering and reaction proceed through the same compound state of total 
momentum, J', and are dominated by one incident partial wave, 2’. 
With this assumption, and substituting I =0, i = 1/2, s' =1/2 in 


Eq. (12) we obtain 


" t 
2 Bes Ba ; 
a 1 e ' ' 1 t a = ~ 
AB (a3; 3) =-y 2(2 RAE a6 cos IL) 7 Or g't (BS) 
sin 2 Opin 
2 
where J' must be either £' +4 or 4£' - 3, depending on the nature 


of the compound resonance. By substituting Eq. (B5) in Eqs. (9) to (11), 


one obtains 
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where is the measured reaction cross section. 


ae 
(cy. Se 
The optical model indicates a maximum in the s wave strength 


19 


18 
function at A*~ 56, and the analysis of Tucker, et al., ~ shows 


that s wave dominance of the outgoing neutrons is consistent with the 


6 
? Knowing that the excited 


measured inelastic cross section in Fe 
state spin is ay we conclude £' =2 and J' is either 5/2 or 


3/2. For J' = 3/2; 


(o 3/2 ) fel aha 53° 
ee * rf'=2’exp 2 2 
= - 3/2 [1 + re (cos @) Jan (B7a) 
a sin 2ba1 
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For J' = 5/2: 


é 5/e Dae 524? 
A do sot =2bexp ; 1 + Sp (cos 6) + 6 p (cos 0) 
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Using Eq. (16), we obtain the expected cusp effects in the total cross 


section. For J' = 3/2: 


cos 203? 
fea = 1 } (B8a ) 
at fae) ar - sin 2834 3/2 
For J' = 5/2 
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Since we have assumed isolated resonance behavior in the cross sections 
in this case, we can write the phases 6 as explicit functions of energy. 
Defining ¢€51 = (EY 0 AGO 2); where E. is the resonance energy 


and Dos is the full width at half maximum, we find 


! 
2 sin” 57 = H+ é) 


1 


cos 2 57: == (1 = aa + eo, (BO ) 
~ ginems- eae e/(1 + >) ; 
where € ranges from - © to +o across the resonance. 
Below-threshold effects are possible in both the differential 


elastic and total cross sections if the special assumptions are valid. 


Unfortunately, this cannot be checked in the present experiment, since 





120 


another resonance feature appears just below threshold, whose spin is 
unknown. 

The above-threshold feature can be analyzed for possible cusp 
effects within the framework of these assumptions, by using Eqs. (B7) 


with Eqs. (B9) substituted. 


es) ( 2 } 
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(B10) 
(1 + P,,(cos 0), if J' = 3/2 
«| 


js 
(1 + (8/7) P(cos @) + (6/7) P)(cos 6)), if J’ = 5/2 


a ees 


ATS egy ey ee 1 +) 


(B11) 


Equations (B10) and (B11) show an interesting property. The 
cusp term in expressions for elastic cross sections (differential and 
integrated) is negative over the entire resonance. The total cross 
section expression is negative for - l<e¢ée<1, and positive for 
e<-l1 ands+eo>dl. This can be interpreted simply as an increase 
in the total cross section resonance width due to the reaction. How- 
ever, the experimental neutron energy spread (~ 7 kev) was too broad 
to allow measurement of such a difference in the widths of the elastic 
and total cross section resonance features. 

In conclusion, we remark that the below-threshold feature (see 
Fig. 14, for example) which we have described as ‘mirrored’ across 
threshold is just the kind of cusp effect one might expect from the 
theory. We were disappointed that our cusp analysis did not demonstrate 


this feature to be a cusp anomaly. 


APPENDIX C 
A PROOF 


We wish to show that for the case of many resonances in the 


averaging interval one finds under certain conditions 
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where 6 is the Kronecker delta. Expanding the left side of (Cl), 
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lize 


Then expression (C2) becomes 
= ( Aa + Bb) (C3) 


We now assume, for £#2' and/or J #J' and within an averaging 
interval, that the quantities A, B, a, and b may be written as the 


sum of a mean value and a fluctuating part; for example, A =A° +AA, 


where ( AA )=0O. Thus, expression (C3) may be written 
( p°a® + Bd? + A°A.a + a°AAt+ BAKDTDAB+AAAARAt+ABAD) (Ch) 


Meyerhof showed in Ref. / that 


1 


y? - cos 28 
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(B)#0O (20) 


Also 
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cos 28 
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on the same basis as the proof of (20). Therefore, A° =B° ~O and 


expression (C4) is reduced to 
(AAAatABAb) (C5) 
where we have also used ( a®°AA) = ( 5a ) = 0. 


1 
QA and AB are functions of 5 : While Aa and Ab refer only 
J i 
to 5): Hor of f L and/or J # J , the fluctuation over the averag- 


ing interval in functions of 5 


J are uncorrelated to those in func- 


t 
tions of 8)1; thus, expression (C5) is equal to zero, and we have 


proved Eq. (C1) for this case. 
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For J=J' and £2 =2' the IMS of Eq. (Cl) reduces to 


j 
a ane ue 
J L 
(g 4 iam 
re nee 

sin 2 O 1 
and we may use the result proved in Ref. 7, 


oC Le 52 


Je 2 (on) 
“it Lan 2 of ine 0 j (9) 


which completes the proof of Eq. (Cl). 


APPENDIX D 
GAMMA RAY SUPPRESSION PREAMPLIFIER 


A space charge saturation gamma-ray suppression preamplifier 
was used to distinguish neutron induced from gamma-ray induced pulses 
in the detector. The circuit is shown in Fig. e2. It is basically 
that of re, modified to include more external controls, a diode 
clipper on the nonlinear output, and an integrating network for the 


remaining positive overshoot. The external control on resistor R 


53 


1 
was suggested by the work of Funsten and Cobb. 


Certain organic scintillators emit light under particle bom- 
bardment, which is composed of an initial spike of fluorescence with 
about 50mu sec decay time, followed by one or more slower components 
of greater than 200 mp sec decay cc Ion recombination may ac- 


oe, 


count for the slow components. The ratio of the amplitude of the 
initial spike to that of the first slow component is a function of the 
type of particle causing the fluorescence. 

The circuit in Fig. 22 is designed to produce linear voltage 
pulses, 1.e., voltage pulses which are proportional to the amplitude 
of the initial fluorescence, at dynode. 10. The remainder of the re- 
sistor chain is arranged deliberately to produce space-charge satura- 
tion in the region between dynode 14 and the anode, during the initial 


spike. When the slow-component cascade electrons arrive in the region 


of space-charge saturation they are sufficiently repelled to produce a 
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Bitgemce . 


Circuit diagram of a space charge saturation gamma 


ray suppression preamplifier. The four externally controlled elements 


are marked "FOC, 


" "acc," "R,," and "DYN 14." "LIN" means linear 


pulse output. "y supp means the gamma ray suppression (nonlinear) 


output. 
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positive overshoot pulse at dynode 14. Thus, the height of the positive 
overshoot is a measure of the ratio of slow component to initial spike 
amplitude and it is possible to use this effect to distinguish the type 
of particle which initiated the fluorescence. 

One is particularly interested in distinguishing neutron- 
initiated pulses (i.e., from recoil protons) from pulses caused by 
gamma-rays (i.e., due to electrons), and the method has proved useful 
for this purpose. Our particular circuit is arranged such that the 
negative pulse at dynode 14, caused by the initial spike for either in- 
cident neutrons or gamma-rays, is clipped by a diode. The remaining 
positive overshoot is, of course, no longer a linear effect, and we 
have chosen to integrate this pulse in an attempt to maximize the dis- 
tinction between neutron and gamma-ray initiated pulses. 

The circuit was effective in gamma suppression for neutron 
energies greater than about 0.3 Mev; however, some care must be taken 
in biassing out gamma-ray induced pulses and to select neutron caused 
pulses, particularly when working at low neutron energies. 

An example of the nonlinear spectrum obtained at 1.7-Mev neu- 
tron energy is given in Fig. 23. A nonlinear spectrum of gamma-ray 
initiated pulses is also shown for comparison. These gamma-rays were 
produced under Van de Graaff machine operating conditions with the 
proton beam incident upon a stopper in front of the tritium target. 

While working at low neutron energies (<0.5 Mev) we found that 


very small dynode 14 voltages were necessary to provide good gamma-ray 
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Fig. 23. Nonlinear pulse spectra obtained at E =e 7 Men 
from a space charge saturation gamma ray suppression preamplifier. 
The detector crystal was in the 0° neutron beam direction. The 
solid circles show the nonlinear spectrum with the proton beam on 

3 3 | 
the T° target (T(p,n)He~). Open circles show the nonlinear spectrum 
with the proton beam blocked by a tantalum beam stopper in front of 


the T° target. 
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suppression. This produced smaller nonlinear output pulses, which 
could be compensated by increasing the high voltage on the photomulti- 
plier, but which eventually represented an essential limitation on the 


circuit function because of the increasing influence of tube noise. 





APPENDIX E 
MULTIPLE SCATTERING AND FINITE SCATTERER CORRECTIONS 


The slab scatterers used in the experiment were between 1/5 
and 1/2 neutron mean free paths in thickness. Thus a significant 
number of multiply scattered neutrons were detected. Cross section 
data must be corrected for this effect. A second correction was re- 
quired because the detection crystals were close enough to the slab 
scatterer to detect once-scattered neutrons over an angular range of 
about +8 degrees in polar angle with respect to the line of crystal-- 
scatterer centers. One can see this by examining Fig. 1 which shows 
the scatterer and detectors in typical positions. Because the detectors 
were sampling the differential cross section over this angular spread, 
peaks in the experimental angular distribution were reduced and valleys 
filled. We will refer to the correction for this effect as a ‘finite 
scatterer correction.” 

We were not prepared to make the computer analysis necessary 
for obtaining accurate corrections because for cusp analysis a first 
order estimate was considered satisfactory. 

The measurements of differential cross sections versus energy 
were made for only one or two angles at a time. Hence, the experimental 
points were correlated in energy but not in angle. The absolute error 
in the (uncorrected) measured cross section was approximately +15 per- 


cent (see Error Analysis) and was independently established for each 


lee ad 
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angle. Therefore, when one cross plots to present the angular distri- 
bution for different energies, the points determining the shape of the 
angular distribution are separately uncertain by #15 percent. 

The shape of the differential scattering cross section is, of 
course, essential to a determination of the multiple scattering correc- 
tion. We used first order multiple scattering and finite scatterer 
corrections in this experiment in two ways. We corrected the optical 


ml for the angular distribution of cerium 


model prediction of Moldauer 
and compared the results with experimental data (see Fig. 8). We also 
corrected the calculated "“energy-averaged" cusp effects before applying 
them to the experimental data for Ce, Ba, and Zr. 

In the Fe experiments, the structure in the data did not justify 
any detailed estimate of the corrections, short of a full computer study 
on the angular and energy dependence of multiple scattering. On the 
basis of characteristic scatterer dimensions we estimated an upper 
limit of 2O percent for the corrections in the case of iron. 

In the lithium experiments, the scatterer was 0.17 neutron 
mean free paths in characteristic dimension. This dimension is small 
in comparison with other scatterer thicknesses typically used in this 
experiment. The angular distribution is peaked backwards and the lab 
neutron energy loss for back scattering is from 30 to 4O percent of 
incident energy (below the detection bias of the experiment). All of 


these factors, particularly the last, contribute to an estimate that 


the correction for multiple scattering in this case is small. The 
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general agreement in magnitude between the measured lithium angular 
distribution and other work is shown in Fig. 9 and supports this con- 
clusion. 

It remains to discuss in detail the first order multiple scat- 
tering and finite scatterer corrections for Ce, Ba, and Zr. We used 


36 


a procedure and notation presented by Nauta. Three assumptions 


underlie the calculation: first, the ratio is independent 


Si 41/Sq 
of K, the order of scattering, where ok is the probability of scat- 
tering elastically K times before leaving the scatterer; second, the 
angular distribution of multiply scattered neutrons is the same as in 
an infinite medium of the same material; and third, there is no loss 
of lab energy upon scattering. The scatterer is characterized by a 
parameter & which is the fraction of neutrons which, having been once 
scattered elastically, have undergone at least one more interaction. 
Neglecting third and higher order scattering, the necessary 
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formulas are: 


Corr fg (9) do. = [1+ (S,, (6)/s, (e))] do. (9) (E1) 

where 
i el 7 2 
s.°"(6)/8,°"(@) = (2x 8/o,,do,_) £ [2/(2i +1))(a,)~ P.(cos 6) 
e dl ower : i ai 
=O 
1 il. 
a, = Bed + 1) J do.7 P, (u)dy 

and where ADie is the true ditfermentialseress sections. sdo is 


exp 


the measured differential cross section; o,, is the total cross section; 


dé 
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el el ; el el 
and do _jaeae, (S, (@))/s, yo Withee SS = f 55 (@)aQ. The last 
definition emphasizes the fact that these equations are circularly 
related and that approximations of some kind are necessary to obtain 
a solution. Our procedure was to substitute an optical model angular 
aff 


distribution for do. : 


v thus we calculated a correction to the opti- 


cal model prediction for comparison with experiment. Satisfactory .. 
agreement between the corrected optical model distributions and experi- 
ment was considered justification for applying these corrections to 

the comparison of predicted cusp effects with experiment. This general 
procedure is analogous to that usually adopted in solving Eq. (El); 


one assumes a function dg,, and by iteration of Eq. (E1) forces it 


di 


to converge to nS gees We stopped this process at a first iteration. 

In essence, the method separates that part of the multiple scat- 
tering which depends upon the angular distribution characteristic of 
the scattering material from a part which depends upon the size, shape, 
and aspect (with respect to the beam axis) of the scatterer. The 
latter dependence is contained in the parameter 5. 

Lacking a better estimate of 8, we assumed that the quantity 
?,/?1, in the detailed work of A. Langsdorf, Jr., R. O. Lane, and 
caer, . Menahan- ' is 1-8, where ?) is the emergent singly scattered 
flux per unit solid angle and P11 is the unit idealized scattered 
flux per unit solid angle. The parameter & is dependent upon the 
angle of scattering at which the cross section is measured because 


36 


slab scatterers are extremely asymmetric. In the work of Nauta, who 


used ring scatterers, 68 could be assumed independent cf angle. 
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Accordingly, we used the values of 9,/91, piven an rae. 9 of 
Ref. 37 to obtain a §& for each scatterer and for each angle of inter- 
est. In most cases extrapolation of the curves was necessary. These 
curves are entered with a particular value of h, the scatterer thick- 
ness, in units of neutron mean free path. 

Our scatterer thicknesses, calculated at the neutron energies 
of interest, were: for Zr, H = 0.32; for Ce, h = 0.46; and for Ba, 
|e Get OS we 

The moments a. of the optical model angular distribution were 
included to ay: Multiple scattering corrections as a function of 
angle were then computed using Eq. (El) and the appropriate values of 
5. The degree to which corrected optical model predictions agreed 
with experimental data is indicated in Fig. 8. We estimated the over- 
all uncertainty in the absolute multiple scattering corrections result- 
ing from the assumptions and approximations in the foregoing treatment 
to be £30 percent. Expressed as the fraction en which is con- 
venient for modifying data, the uncertainty in multiple scattering 
corrections adds to the uncertainty obtained in the experimental error 
analysis. The uncertainty in fs. varies with scatterer and detec- 
tion angle, as indicated in Table III. 

We now consider the calculation of finite scatterer corrections. 
At each mean polar angle, 64 once scattered neutrons are detected 
over an angle spread of ae + AO. The magnitude of 2 AO was de- 
termined graphically from the geometry of the experiment, see Fig. l, 


and was found to vary from 10 to 35 degrees. The angular spread is a 
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funetion,of @. because the scatterer-to-detector distance varied 
depending upon the experimental situation (counting rate, background to 
effect ratio) and also because the projected area of the scatterer, 
normal to the scatterer-detector line of centers, varies with oe 

As in the correction for multiple scattering, the procedure 
was to correct an optical model prediction for experimental effects. 


The finite scatterer correction is defined by 


ao worr ~ fois .69) ad xp (E2) 


— (@) was computed using 


N 
fee) me do. /(1/N Sat do, .4(sin @_/sin 6) (E3) 


where do.4 69) is the optical model prediction at mean angle 6, 
and the denominator represents an average of the optical model distri- 
bution over N equal divisions of the spread 2 AO weighted by the 
solid angle subtended by the nth angle division. N was selected to 
give a representative average and varied from 5 to 8 depending upon 
the size of @2 a GC. 

The results of this calculation are given in Table III. The 
magnitude of typical finite scatterer corrections is indicated on the 
angular distribution shown in Fig. 8, where it can be seen that the 


correction is small except in the valley near a = o0>. 
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TABLE III. Multiple Scattering and Finite Scatterer Corrections 


Uncertainty 
@ fi fa in f pb 
O m.s. TAS. m.S. rss 
(percent ) 
Cerium 30° 0.95 1.04 
50° Gms 0.98 
O° 0.80 0.92 
90° 0.75 0.90 10 
iDlerg 0.78 0.99 9 
130° 0.80 1 ere 7% 
15 Oy 0.85 1.05 6 
Barium con 0.65 0.90 18 
125° OA 70 OO 14 
Zirconium 60° 0.90 ~1.0 4 
130° 0.85 ~] 6 
Tron MBE OS 10 ~1.0 ~10 
Lithium Adal ~1.0 1-0 ~ 5 





The factors in this table may be applied to the::experimental 


cross section data presented in this report, by using 


do = of £ do ‘ 
Corr M.S.) “ies. exp 
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